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Abstract  .  ~  - . ' 

^  In  this  report,  we  give  a  brief  summary  of  the  work  done  py  the  Pis  and  their 
students  on  the  design  and  analysis  of  robust  feedback  systems/  Under  this  contract, 
we  have  conducted  research  on  a  number  of  problems  in  system  and  control  theory. 
In  particular,  we  have  obtained  significant  new  results  on  nj,  control  theory,  robust 
stability  analysis  of  feedback  systems,  robust  stabilization  of  uncertain  systems,  ap¬ 
proximation  and  control  of  infinite  dimensional  systems,  etc.  In  this  report,  we  will 
present  a  brief  summary  of  the  main  results  obtained  under  the  AFOSR  contract.  / 


*. 


1  Introduction 

Our  (Pis’  and  their  students’)  research  work  in  systems  and  control  has  been  supported 
for  the  last  3  years  by  the  AFOSR.  We  have  made  significant  contributions  to  H#,  control 
theory,  robust  control  -  analysis  and  synthesis,  control  of  infinite-dimensional  systems,  control 
of  time-varying  systems,  etc.  Our  work  on  these  problems  naturally  divides  up  into  2  parts: 
work  done  Professor  Khargonekar  and  his  group  and  the  work  done  by  Professor  Tannenbaum 
and  his  group.  While  they  have  addressed  related  problems,  their  work  and  approaches  are 
complement  each  other,  and  it  is  easiest  to  describe  the  research  accomplishments  under  the 
AFOSR  contract  into  two  separate  parts. 

We  have  published  most  of  the  work  under  the  contract  in  various  journals  and  conferences 
as  can  be  seen  in  the  attached  list  of  publications.  For  reasons  of  space  and  time,  we  will 
only  present  brief  descriptions  of  the  results  and  give  appropriate  references  where  the  results 
may  be  found.  We  have  strictly  restricted  the  discussion  to  our  own  work  and  references  to 
our  own  papers  since  a  complete  discussion  of  related  references  will  make  this  report  too 
long  and  unwieldy.  The  interested  reader  is  referred  to  our  published  papers  where  extensive 
discussions  and  lists  of  related  references  may  be  found. 


2  Work  done  by  Professor  Khargonekar  and  his  group 

2.1  State-Space  Approach  to  Hoo  Control  Theory 

Motivated  by  robustness  considerations,  Zames  introduced  the  problem  of  #<*>  optimal  con¬ 
trol.  The  essential  idea  was  to  design  a  controller  to  optimize  performance  for  the  worst 
exogenous  input.  Thus,  while  in  Kalman  filtering  and  the  LQG  problem,  the  power  spec¬ 
trum  of  the  exogenous  input  (noise)  is  assumed  known  (usually  white  noise),  in  the 
control  problem  the  power  spectrum  is  assumed  unknown  and  the  controller  is  designed  for 
the  worst  case. 

Early  research  in  Hoo  control  theory  was  conducted  using  frequency  domain  methods. 
The  key  tools  were  the  Youla-Jabr-Bongiorno-Kucera  parametrization  of  all  stabilizing  con¬ 
trollers,  inner-outer  factorizations  of  transfer  functions,  Nevanlinna-Pick  interpolation  the¬ 
ory,  Nehari  distance  theorem,  the  commutant  lifting  theorem,  etc.  This  frequency-domain 
operator  theoretic  approach  is  currently  a  very  active  area  of  research  and  a  number  of  dif¬ 
ferent  avenues  are  being  explored.  The  work  done  by  Professor  Tannenbaum  and  his  group 
is  along  these  lines  and  a  summary  of  results  obtained  along  this  approach  may  be  found 
elsewhere  in  this  report. 

A  major  new  development  in  H^  control  theory  in  the  last  two  years  has  been  the 
introduction  of  state-space  methods.  This  has  led  to  a  rather  transparent  solution  to  the 
standard  problem  of  H9 0  control  theory.  This  solution  is  remarkably  similar  to  the  classical 
LQG  solution  with  appropriate  differences  that  reflect  the  differences  in  the  Hoo  and  the 
LQG  performance  criteria.  As  a  result  it  is  possible  to  give  a  simple  interpretation  to  the 
so- called  Hoo  central  controller.  Analogous  to  the  LQG  theory,  the  solutions  are  given  in 
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terms  of  solutions  to  algebraic  and  differential  matrix  Riccati  equations. 

Consider  the  finite-dimensional  linear  time-invariant  system  E 

dx 

—  =  Fx  +  G\tv  +  G2u, 

at 

z  —  H\X  -j-  Ji\iv  +  J\ 2U, 

y  =  H^x  +  J21W  +  J22U  (1) 

Here  x,  w,  u,  z,  y  denote  respectively  the  state,  the  exogenous  input,  the  control  input,  the 
regulated  output,  and  the  measured  output.  Kalman  made  pioneering  contributions  to  the 
problem  of  designing  a  controller  for  optimizing  the  variance  of  z  when  w  is  a  stochastic 
process.  This  classical  solution,  commonly  known  as  the  LQG  theory,  has  had  a  very 
significant  influence  on  linear  multivariable  control  theory  over  the  last  30  years. 

The  standard  problem  of  Hoo  control  theory  is:  Given  the  linear  system  E  and  a  positive 
number  7,  find  a  causal  (dynamic)  controller  K  such  that  the  closed  loop  system  is  well- 
posed,  internally  stable,  the  closed  loop  input-output  operator 

Txw:  L2  — ♦  L2:  w  z 

is  bounded,  and  the  induced  norm 

\\TZW\\  :=  suP(jl~j^  :  IMI2  7*  0}  <  7- 

A  controller  K  is  called  admissible  iff  it  is  causal,  the  closed  loop  system  is  well-posed 
and  internally  stable.  If  K  is  also  an  FDLTI  system 

§  =  M  +  By, 

u  =  Ci  +  Dy  (2) 

then  the  closed  loop  system  it  we’J-posed  iff  (/  —  J22D)  is  invertible  and  the  closed  loop 
system  is  internally  stable  iff  the  sed  loop  system  matrix 

'  F  +  G2(I  -  DJ22)~1DJ22  G2(I  -  DJ22)~'C 
B{I  -  J22D)-lH2  A  +  BJ22{I  -  DJ22)~lC 

has  no  eigenvalues  in  the  closed  right  half  complex  plane.  Moreover,  in  this  case,  if  we 
let  Txw(s)  denote  the  closed  loop  transfer  function  matrix,  then 

l|r~l  =  »rre|L  :=  sup{<7(TU*))  ■■  Se(»)  >  0}. 

It  is  a  well  known  fact  that  under  suitable  assumptions,  the  LQG  controller  is  also  the 
unique  solution  to  the  problem  of  minimizing  the  quadratic  norm  of  the  closed  loop  transfer 
function  matrix  Txw  over  all  internally  stabilizing  controllers.  Thus,  the  key  difference  be¬ 
tween  the  LQG  and  the  H0 0  control  problems  is  in  the  choice  of  norm  on  transfer  functions. 
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This  is  intimately  related  to  the  underlying  assumptions  on  the  exogenous  signals  as  men¬ 
tioned  in  the  introduction.  Also,  it  turns  out  that  the  H^  norm  axises  naturally  in  many 
robust  control  problems. 

In  order  to  present  the  most  important  concepts  clearly,  we  will  make  certain  simplifying 
assumptions.  Most  of  these  assumptions  can  be  easily  removed. 

A. I  Jn  =  Jii  —  0* 

A.2  f„  [  Ht  Jn  ]  =  [  0  /  ]  . 

A.3  Jn  [  C,  fu  ]  =  [  0  /  J  . 

A.4  (F,  G\,H\)  is  stabilizable  and  detectable. 

A. 5  (F,  G2,  Hi)  is  stabilizable  and  detectable. 

A.6  The  controller  K  is  an  FDLTI  system. 

Assumption  A.l  implies  that  there  is  no  direct  transmission  from  the  exogenous  input  to 
regulated  output  and  from  the  control  input  to  the  measured  output.  The  latter  assumption 
ensures  that  any  proper  rational  matrix  K(s )  leads  to  a  well- posed  feedback  system.  As¬ 
sumption  A.2  is  quite  common  in  the  LQG  literature  and  amounts  to  assuming  that  there  is 
no  cross  term  in  the  formula  for  ||z||2  and  the  penalty  on  the  control  input  u  is  normalized, 
i.e., 

\\z\\2  =  x'H'1H1x  +  u'u. 

Assumption  A.3  is  the  dual  of  assumption  A.2  and  is  analogous  to  the  standard  assump¬ 
tion  in  the  Kalman  filtering  problem  that  the  process  noise  and  the  measurement  noise  are 
uncorrelated  and  that  the  measurement  noise  is  nonsingular  and  normalized  .  Assumption 
A.4  is  a  technical  assumption  and  guarantees  that  the  corresponding  controllers  are  admis¬ 
sible.  Assumption  A.5  is  necessary  and  sufficient  to  guarantee  the  existence  of  an  internally 
stabilizing  controller  for  the  system  E. 

The  assumption  A.6  also  causes  no  loss  of  generality  as  shown  by  Khargonekar  and  Poolla 
[1989].  They  showed  that  the  infimum  of  the  norm  of  Tzw  over  all  causal  internally  stabilizing 
nonlinear  controllers  is  no  less  than  that  over  all  FDLTI  internally  stabilizing  controllers. 

2.2  The  State>Feedback  H Control  Problem 

We  will  first  consider  the  state-feedback  Hoo  problem.  In  other  words,  in  the  system  E,  let 
us  assume  that 

H2  —  I  and  J2 1  =  0.  (3) 
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(It  should  be  noted  that  the  state-feedback  problem  does  not  satisfy  the  assumption  A. 3 
above.)  Under  the  assumption  Jn  =  0,  Petersen  considered  the  problem  of  finding  a  real 
matrix  L  such  that  F  +  GiL  is  asymptotically  stable  and 

\\Hi(sI  —  F  —  GjZ/)~1G1||00  <  7. 

He  obtained  a  necessary  and  sufficient  condition  for  the  existence  of  L  in  terms  of  an 
algebraic  Riccati  equation  with  an  indefinite  quadratic  term.  This  was  quite  an  interesting 
result  and  led  to  a  number  of  subsequent  developments. 

In  view  of  the  fact  that  frequency-domain  control  theory  results  usually  led  to  (ap¬ 
parently)  high  dimensional  controllers,  Khargonekar,  Petersen,  and  Rotea  [1988]  investigated 
whether  static  gains  were  actually  optimal.  This  problem  was  motivated  by  the  one  of 
Kalman’s  key  insights  that  control  is  a  function  of  the  state.  Thus,  we  expected  the  following 
result  to  be  true: 

Theorem  2.1  Consider  the  system  E  and  suppose  (3)  holds.  Then 

7*  =  inf  {  ||7hJoo  :  K  is  an  internally  stabilizing  dynamic  state-feedback  controller  } 

=  inf  {  ||7’IU)||l0  :  K  is  an  internally  stabilizing  static  state-feedback  gain.  } 

Theorem  (2.1)  was  proved  by  Khargonekar,  Petersen,  and  Rotea  [1988]  under  the  con¬ 
dition  that  Ju  =  0.  It  was  generalized  further  by  Khargonekar,  Petersen,  and  Zhou  [1987], 
Zhou  and  Khargonekar  [1988],  Doyle,  Glover,  Khargonekar,  and  Francis  [1989]. 

The  following  theorem  shows  how  one  can  obtain  state-feedback  control  laws  by  solv¬ 
ing  certain  algebraic  Riccati  equations.  It  is  taken  from  Doyle,  Glover,  Khargonekar,  and 
Francis  [1989].  Results  analogous  to  Theorem  (2.2)  were  also  obtained  in  our  earlier  papers 
Khargonekar,  Petersen,  and  Zhou  [1987],  Zhou  and  Khargonekar  [1988]. 

Theorem  2.2  Consider  the  system  E.  Suppose  assumptions  A.l,  A. 2,  and  (3)  hold.  Then 
there  exists  an  admissible  controller  K  such  that 

Halloo  <  7  (4) 

if  and  only  if  there  exists  a  (unique)  symmetric  matrix  P  such  that 

F’P  +  PF  +  P(£0k-G2G'2)P  +  H’lHi  =  0,  (5) 

F  +  (Gffi  —  GjGj)/3  is  asymptotically  stable,  and  P  >  0.  In  this  case,  the  control  law 

u  =  —G'2Px  (6) 

internally  stabilizes  E  and  (4)  holds. 
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There  are  a  number  of  interesting  features  to  this  result.  Note  that  the  control  law  is 
obtained  by  solving  an  algebraic  Riccati  equation  which  is  analogous  to  the  classical  results  of 
Kalman  on  the  linear-quadratic  optimal  control  problem.  The  algebraic  Riccati  equation  (5) 
is  similar  to  the  corresponding  equation  in  the  linear-quadratic  optimal  control  problem 
except  that  the  quadratic  term  in  (5)  is  indefinite.  Indeed,  (5)  is  identical  to  the  ARE  that 
arises  in  linear  quadratic  optimal  game  problems.  This  is  intuitively  appealing  since  in  the 
Hoc  control  problem  the  inputs  to  and  u  act  as  opposing  players:  the  exogenous  input  to  tries 
to  maximize  (the  norm  of  )  of  z  while  u  is  designed  to  minimize  it.  This  connection  between 
linear-quadratic  games  and  control  theory  has  been  discussed  and  explored  in  a  number 
of  recent  papers.  See,  e.g.,  Khargonekar,  Petersen,  and  Zhou  [1987],  Limebeer,  Anderson, 
Khargonekar,  and  Green  [1989],  and  the  references  therein. 

The  paper  by  Limebeer,  Anderson,  Khargonekar,  and  Green  [1989]  contains  analogues  of 
Theorem  (2.2)  for  linear  time-varying  systems. 

Note  that  as  7  goes  to  00,  the  controller  (6)  approaches  the  LQR  solution.  In  other 
words,  as  the  norm  constraint  on  the  closed  loop  transfer  function  is  relaxed,  the  control 
law  given  by  (6)  approaches  the  corresponding  LQR  control  law. 

The  results  of  Khargonekar,  Petersen,  and  Zhou  [1987],  Zhou  and  Khargonekar  [1988] 
also  apply  to  the  general  singular  case  (J21  is  not  necessarily  full  rank).  In  these  papers 
the  singular  case  is  taken  care  of  by  introducing  certain  small  perturbations  to  make  the 
problem  nonsingular.  Since  we  are  dealing  with  a  strict  inequality  in  (4),  the  existence  of 
an  appropriate  small  perturbation  is  not  difficult  to  establish. 

Theorem  (2.2)  naturally  leads  to  the  Hoo  analog  of  the  inverse  problem  of  optimal  control 
formulated  by  Kalman.  This  inverse  problem  of  Hoo  control  theory  was  considered  by  Fujii 
and  Khargonekar  [1988]  who  showed  that  state-feedback  //«,  controllers  inherit  the  nice 
robustness  properties  of  the  LQR  controllers  and  in  a  certain  sense  are  even  more  robust. 

2.3  The  Hoo  Filtering  Problem 

In  the  LQG  theory,  state-feedback  linear-quadratic  optimal  control  and  Kalman  filtering 
problems  turn  out  to  have  dual  structures.  The  situation  is  somewhat  similar  in  the  Hoo 
case.  Filtering  problems  were  first  considered  by  Doyle,  Glover,  Khargonekar,  and  Francis 
[1989]  under  the  assumptions  that  (i)  the  process  is  stable,  and  (ii)  the  initial  state  in  known. 
Recently,  Nagpal  and  Khargonekar  [1989]  have  developed  a  fairly  complete  estimation 
theory.  Here  we  will  give  a  brief  summary  of  some  of  their  results. 

Consider  the  FDLTI  system 

I  =  fl+G-w' 

z  =  H\x, 

y  =  H2x  +  J2iw.  (7) 

In  this  section,  we  will  assume  that  (F,  G\)  is  reachable,  (F,  H2)  is  detectable,  and  A. 3  holds. 

The  problem  is  to  estimate  the  output  z  using  the  measurements  y.  Nagpal  and  Khar¬ 
gonekar  [1989]  have  considered  both  the  filtering  (estimator  is  required  to  be  causal)  and  the 


smoothing  (estimator  is  allowed  to  be  noncausal)  problems.  So,  let  T  be  a  causal  estimator 
and  let 


*  =  Hv)- 

The  estimator  T  is  called  unbiased  if 
y(r)  =  0  Vt  <  t  =»  z(t)  =0  Vr  <  t. 


Let  us  now  define  cm  H0 0  i.e.,  a  worst  case,  performance  measure.  There  are  two  cases  to 
consider:  (i)  initial  state  is  known  (and,  without  loss  of  generality  =0),  and  (ii)  the  initial 
state  is  unknown.  In  the  latter  case,  it  is  assumed  that  the  best  a  priori  estimate  of  x(0)  is 
zero.  In  other  words,  the  estimator  initial  state  is  taken  to  be  zero.  Let  0  <  T  <  oo.  Define 


z  —  z 


J\{F,  T)  :=  sup{^—~ M  :  to  €  L2[0,T],0  #  ||u,||2,x(0)  =  0}, 
w  2 


(8) 


and 


z  -  z 


J2(JrfB,T)  :=  sup{ -  4'*—  r^-n/2  :  w  e  T2[0,T],x(0)  =  x0,  \\w\\22  +  x'0Rx0  ^  0}(9) 

|«^||2  +  X0J?Xo]1/2 


In  the  definition  of  J2,  R  is  taken  to  be  positive  definite.  The  key  difference  between  J\ 
and  J2  is  that  J2  mesures  the  worst  case  performance  over  all  possible  w  and  x(0).  The 
expression  [||u>||2  +  XqRx0]1/2  is  a  total  measure  of  energy  in  the  exogenuos  variables  w,  x(0). 
The  following  results  are  taken  from  Nagpal  and  Khargonekar  [1989]. 


Theorem  2.3  Consider  the  FDLTI  system  (7).  Let  7  >  0  be  given  and  let  T  =  00.  Then 
there  exists  an  unbiased  linear  filter  such  that 

J2(F,R,T)  <  7  (10) 

if  and  only  if  there  exists  a  bounded  symmetric  matrix  function  Q(t),  t  E  [0, 00)  such  that 

9(f)  =  FQ{t)+Q(t)F'  -  0(0)  =  JT',(11) 

and  the  system 

4(0  =  [/•  -  -  — ~^-))e(i) 

7 

is  exponentially  stable. 

Moreover,  if  these  conditions  hold  then  one  filter  that  satisfies  the  performance  bound 
(10)  is  given  by 

i(()  =  Fi(t)  +  Q(l)H'Mt)  -  i(0)  =  0,  (12) 


i(()  = 


(13) 
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Theorem  2.4  Consider  the  FDLTI  system  (7).  Let  z(0)  =  0.  Let  7  >  0  be  given  and  let 
T  =  00.  Then  there  exists  on  unbiased  linear  filter  such  that 

J1{Jr,T)<i  (14) 

if  and  only  if  there  exists  a  real  symmetric  matrix  Q  such  that 

FQ  +  QF'  -  QH2H2Q  +  9*L!h9.  +  G&  =  0,  (15) 

7 

the  matrix  (F  —  Q(H2H2 - b-5-))  has  no  eigenvalues  in  the  closed  right  half  plane,  and 

Q>  0. 

Moreover,  if  these  conditions  hold  then  with  Q(t)  replaced  by  Q,  the  filter  given  by  (12) 
and  (13)  satisfies  the  performance  bound  (14)- 

These  results  are  the  H, »  analog  of  the  well  known  Kalman  filtering  results.  The  key 
differences  are  in  the  Riccati  differential  equation  (11)  and  the  algebraic  Riccati  equation 
(15)  which  are  very  similar  to  the  covariance  equations  for  the  Kalman  filtering  problem  with 
the  exception  of  the  term  QH[H\Q/‘y 2.  Thus,  in  Hoo  filtering,  the  states  to  be  estimated 
influence  the  filter  itself  unlike  in  Kalman  filtering  where  the  optimal  estimate  of  any  state- 
functional  is  obtained  from  the  optimal  state-estimator.  Note  that  as  7  — ►  00,  the  H^  filters 
approach  the  standard  Kalman  filters. 

There  is  another  qualitative  interpretation  of  the  term  QH[H\Qfa2.  It  can  be  regarded 
as  additional  process  noise.  As  a  result,  the  H°°  filter  has  robustness  (in  the  sense  of 
satisfying  (10)  or  (14))  to  variation  in  the  spectra  of  the  exogenous  signals  at  the  expense 
of  performance  if  the  exogenous  input  w  is  actually  a  zero  mean  Gaussian  white  stochastic 
process. 

The  paper  by  Nagpal  and  Khargonekar  [1989]  contains  similar  results  on  Hoo  filtering 
and  smoothing  for  infinite  as  well  as  finite  horizon  cases  for  linear  time-invariant  as  well  as 
linear  time-varying  systems. 

2.4  The  Output-Feedback  H Control  Problem 

Let  us  now  consider  the  output  feedback  problem.  A  key  result  from  the  linear-quadratic- 
Gaussian  control  theory  is  that  in  the  output  feedback  case,  the  optimal  controller  is  the 
Kalman  filter  for  the  optimal  state-feedback  law.  The  above  mentioned  results  on  the  state- 
feedback  Hoo  control  problem  indicate  that  something  similar  may  be  true  in  the  output 
feedback  Ha 0  problem.  The  following  result  is  taken  from  Doyle,  Glover,  Khargonekar,  and 
Francis  [1989]. 

Theorem  2.5  Consider  the  linear  system  £.  Suppose  assumptions  A.1-A.6  hold.  Then 
there  exists  an  admissible  output  feedback  controller  K  such  that  with  the  control  law  u  =  Ky, 
the  closed  loop  transfer  function  satisfies 

Halloo  <  7  (16) 
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if  and  only  if  the  following  conditions  hold: 

1.  there  exists  a  (unique)  real  symmetric  matrix  P  such  that 

C  c' 

F'P  +  PF+  -  G2G'2)P  +  H'1Hl  =  0, 

F  +  (^r1-  —  G2C2)P  is  asymptotically  stable,  and  P  >  0; 

2.  there  exists  a  (unique)  real  symmetric  matrix  Q  such  that 

FQ  +  QF'  +  Q(^-  -  H2H2)Q  +  G\G\  =  0,  (17) 


F  +  Q{^r-  —  H2H2)  is  asymptotically  stable,  and  Q  >  U;  and 
3-  <  1. 

If  these  conditions  hold  then  one  controller  that  satisfies  the  inequality  (16)  is  given  by 

—  =  {F+Q^P]x  +  ZQH'2(y-H2x)  +  G2u 
u  =  - G'2Px ,  (18) 


where 


Z  —  (I 


-.2  1 


(Here  p{.)  denotes  the  spectral  radius.) 

The  notation  in  the  above  theorem  is  meant  to  be  suggestive.  Note  that  the  above  con¬ 
troller  (18)  has  many  similarities  to  (and  some  differences  from)  the  classical  LQG  controller. 
It  is  called  the  central  controller  and  has  many  interesting  interpretations  and  properties. 
As  7  goes  to  oo,  this  controller  approaches  the  LQG  controller. 

It  is  shown  in  Doyle,  Glover,  Khargonekar,  and  Francis  [1989]  that  the  above  controller 
is  an  H<x>  estimator  for  the  state-feedback  control  law  (6)  in  the  presence  of  the  disturbance 

P 

w  =  More  specifically,  let 


r  :=  w  — 


Then  the  equations  for  system  E  become 


dx  G^P, 

—  {F  H - —2 — )x  +  G\r  +  (?2U, 

y  =  H2x  +  J21r. 


(19) 
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Then  if  we  consider  the  problem  of  estimating  the  state-feedback  law  (6) 
u  =  -G2'Px 

and  apply  Theorem  (2.4)  to  the  above  system,  then  the  resulting  filter  is  precisely  the 
output  feedback  controller  in  Theorem  (2.5). 

The  above  result  is  also  closely  connected  to  the  separation  principle  in  the  risk  sensitive 
control  problem.  See  Doyle,  Glover,  Khargonekar,  and  Francis  [1989]  for  further  details  along 
these  lines. 

As  mentioned  previously,  in  Limebeer,  Anderson,  Khargonekar,  and  Green  [1989],  we 
have  extended  these  results  to  linear  time-varying. 

2.5  Robust  Control:  Synthesis  Problems 

The  primary  reason  for  the  use  of  feedback  in  control  systems  is  the  fact  that  almost  all 
control  systems  operate  in  uncertain  environments.  There  are  two  types  of  uncertainties 
that  affect  control  systems:  plant  uncertainty  and  signal  uncertainty.  Plant  uncertainty  arises 
due  to  errors  in  modeling,  changes  in  parameters,  inexact  and  incomplete  data,  modeling 
approximations  while  signal  uncertainty  arises  due  to  exogenous  signals  such  as  disturbances, 
sensor  noise,  etc.  Control  theory  attempts  to  provide  systematic  techniques  and  methods 
to  solve  control  problems  where  uncertainty  is  a  dominant  issue.  Indeed,  from  this  point  of 
view,  robust  linear  control  and  adaptive  control  methodologies  are  different  solutions  to  this 
problem.  Now  one  can  regard  an  adaptive  controller  as  a  nonlinear  time-varying  controller. 
Therefore,  a  fundamental  question  of  interest  is: 

Given  a  feedback  control  problem  for  a  linear  multivariable  system,  what  advantages  (if 
any)  do  nonlinear  time-varying  (NLTV)  controllers  offer  over  linear  time-invariant  (LTI) 
controllers? 

Our  work  with  a  number  of  colleagues  T.  Georgiou,  A.  Pascoal,  K.  Poolla,  R.  Ravi,  M. 
Rotea,  and  A.  Tannenbaum  on  various  aspects  of  this  problem  has  led  us  to  formulate  a 
qualitative  principle  which  captures  known  results  obtained  so  far. 

THE  PLANT  UNCERTAINTY  PRINCIPLE  (Khargonekar-Poolla  [1989]) 

In  robust  multivariable  control  problems,  nonlinear  time-varying  controllers  yield  (signif¬ 
icant)  adavantages  over  linear  time-invariant  controllers  if  and  only  if  there  is  (significant) 
parametric  or  structured  uncertainty. 

The  strong  and  simultaneous  stabilization  problems  have  been  considered  in  many  papers 
in  the  literature.  It  is  well  known  that  in  general,  not  every  LTI  plant  can  be  stabilized  by 
a  stable  controller.  Also,  not  every  pair  of  LTI  plants  can  be  simultaneously  stabilized  by 
an  LTI  controller.  In  contrast  to  the  above  situation,  Khargonekar,  Pascoal  and  Ravi  [1988] 
have  proved  the  following  result: 

Theorem  2.6  Let  F  :=  P\,  Pi, ...,  Pn  be  a  finite  collection  of  finite-dimensional  linear  time- 
varying  (LTV)  plants  such  that  each  Pi;i  =  l,..,n  is  internally  stabilizable.  Then,  F  can 
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be  simultaneously  (internally,  and  in  the  Lp-sense)  stabilized  by  a  stable  finite -  dimensional 
LTV  controller. 

As  a  corollary,  it  was  also  shown  that  any  finite-dimensional  LTI  plant  can  be  internally 
stabilized  by  a  stable,  periodic  controller.  Further,  every  finite  collection  of  finite-dimensional 
LTI  plants  can  be  simultaneously  stabilized  by  a  stable  periodic  controller. 

The  papers  by  Khargonekar  [1989]  and  Khargonekar  and  Poolla  [1989]  contain  a  complete 
discussion  of  known  results  on  the  Plant  Uncertainty  Principle  and  a  number  of  related 
references. 

In  a  parallel  direction,  we  also  investigated  robust  control  synthesis  using  Liapunov  func¬ 
tion  methods  -  ’quadratic  stabilization’  and  robust  control  of  uncertain  systems.  In  Khar¬ 
gonekar,  Petersen,  and  Zhou  [1987,  1990]  we  showed  that  there  is  a  very  deep  connection 
between  i/oo  control  quadratic  stability  and  stabilization  and  Hx  control  theory.  Actually, 
this  connection  goes  back  to  some  early  work  of  Popov,  but  it  had  apparently  been  over¬ 
looked  in  the  more  recent  robust  control  literature.  In  Khargonekar,  Petersen,  and  Zhou 
[1987,  1990],  we  gave  a  complete  and  constructive  solution  of  the  problem  of  quadratic  sta¬ 
bilization  by  LTI  control  by  output  feedback  -  an  unsolved  problem  then  in  the  uncertain 
systems  literature.  Related  technical  results  were  obtained  in  Khargonekar  and  Zhou  [1988]. 
An  interesting  result  supporting  the  plant  uncertainty  principle  was  recently  obtained  by 
Khargonekar  and  Rotea  [1988b],  Rotea  and  Khargonekar  [1989]  where  we  showed  that  for 
systems  with  norm  bounded  uncertainty,  quadratic  stabilizability  via  NLTV  controllers  im¬ 
plies  quadratic  stabilizability  via  LTI  controllers  in  the  state-feedback  case.  Related  results 
were  obtained  in  Rotea  and  Khargonekar  [1988]  where  we  showed  that  if  the  plant  states  are 
available  for  feedback  then  quadratic  stabilizability  by  linear  dynamic  state-feedback  implies 
quadratic  stabilizability  by  nondynamic  state-feedback. 

In  Khargonekar  and  Gu  [1989]  we  have  given  a  new  state-space  solution  to  a  robust 
stabilization  problem  for  systems  with  combined  parametric  and  dynamic  uncertainty  under 
a  minimum  phase  assumption.  This  is  an  important  problem  area  where  very  few  results 
are  currently  available. 

It  is  well  known  that  control  design  invloves  tradeoffs  among  competing  objectives.  In 
Khargonekar  and  Rotea  [1989],  we  have  recently  given  a  complete  solution  to  the  multi¬ 
ple  objective  optimal  control  problem  (in  the  sense  of  Pareto  optimality)  in  the  Hi  norm 
case.  Our  results  include  a  complete  characterization  and  construction  of  all  Pareto  opti¬ 
mal  controllers.  These  are  the  most  general  results  available  on  the  problem  to-date.  The 
computations  invlove  only  solving  algebraic  Riccati  equations  and  cure  easily  implementable 
using  stnadard  matrix  algebra  packages  such  as  MATLAB.  An  interesting  conclusion  from 
this  work  is  that  the  (state)  dimension  of  the  Pareto  optimal  controllers  often  exceed  the 
(state)  dimension  of  the  controlled  plant.  This  is  in  sharp  contrast  to  the  existing  results  in 
the  LQG  and  H control  problems,  where  the  controller  dimension  can  always  be  chosen  to 
not  exceed  that  of  the  controlled  plant.  This  is  a  result  on  the  controllers  complexity  and 
raises  a  host  of  issues  for  future  research. 

Motivated  by  problems  in  adaptive  control,  a  complete  solution  to  a  certain  robust  syn¬ 
thesis  problem  for  positive  real  functions  was  given  in  Anderson,  Dasgupta,  Khargonekar, 
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Kraus,  and  Mansour  [1989].  Roughly  speaking,  given  a  polytope  of  polynomials,  it  is  desired 
to  find  a  fixed  polynomial  (or  a  stable  rational  function)  such  that  the  ratio  of  the  fixed 
polynomial  and  any  polynomial  from  the  given  polytope  is  a  positive  rational  function.  We 
gave  a  necessary  and  sufficient  condition  for  a  solution  to  exist  and  also  a  construction  for  a 
solution  if  it  exists.  This  result  brought  out  an  interesting  connection  to  ihe  growing  body 
of  literature  on  Kharitonov  mathods  for  robust  stability  analysis. 

We  have  obtained  some  key  results  on  the  construction  of  coprime  factorizations  for  linear 
time-varying  systems.  These  are  presented  in  Khargonekar  and  Rotea  [1988].  It  was  shown 
that  a  given  a  linear  time-varying  system  admits  a  coprime  factorization  if  and  only  if  it 
can  be  stabilized  by  dynamic  output  feedback.  In  a  related  direction,  a  generalization  of  the 
Youla  parametrization  for  a  class  of  multirate  digital  control  systems  was  obtained  by  Ravi, 
Khargonekar,  Minto,  and  Nett  [1989].  This  result  may  open  up  anew  research  direction  on 
the  design  of  multirate  control  systems  -  a  very  important  problem  from  a  practical  point  of 
view. 


2.6  Robust  stability  analysis 

In  Barmish  and  Khargonekar  [1988,  1990]  we  formulated  and  solved  a  robust  stability  analy¬ 
sis  problem  for  SISO  systems  containing  uncertain  real  parameters  and  unstructured  uncer¬ 
tainty.  This  paper  still  remains  one  of  a  handful  of  papers  dealing  with  robustness  analysis 
for  simultaneous  parametric  and  unstructured  uncertainty.  The  result  was  an  algorithm 
to  check  the  stability  of  a  SISO  feedback  system  containing  uncertain  real  parameters  and 
multiplicative  unmodeled  dynamic  uncertainty. 

In  a  very  recent  paper  Barmish,  Khargonekar,  Shi,  and  Tempo  [1989],  we  have  discovered 
that  the  robustness  margin  can  be  a  discontinuous  function  of  the  problem  data.  Although 
the  ramifications  of  this  result  will  not  be  clear  for  at  least  a  few  months,  it  certainly  seems 
likely  that  this  will  open  a  new  research  direction  robust  stability  analysis.  In  particular,  it 
raises  the  issue  of  ‘conditioning  properties’  of  robust  stability  margin  problems. 

2.7  Infinite-dimensional  and  sistributed  systems 

In  Gu,  Khargonekar,  and  Lee  [1988,  1989]  we  have  developed  a  new  approach  to  the  problem 
of  finding  finite-dimensional  approximations  for  (possibly  unstable)  infinite-dimensional  sys¬ 
tems  which  converge  in  the  Loo  norm.  Our  approach  uses  the  FFT  algorithm  and  singular 
value  decompositions  and  seems  to  be  computationally  less  demanding  than  the  approaches 
based  on  the  Hankel  operator  theory.  Indeed,  we  have  implented  these  numerical  algorithms 
on  an  IBM  PC,  and  applied  them  to  some  simple  examples,  where  they  seem  to  give  very 
good  results.  They  may  also  have  applications  in  digital  filter  design. 

In  Georgiou  and  Khargonekar  [1989],  we  have  developed  a  new  algorithm  for  spectral 
factorization  of  matrix  valued  functions  using  the  Nevanlinna-Pick  theory.  At  this  point,  the 
computational  potential  of  this  algorithm  is  unexplored  and  it  remains  a  subject  for  future 
investiagtion. 
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1  Introduction 


Referring  to  the  list  of  publications  below  in  Section  7  since  the  start  of 
funding  on  October  1,  1987,  Tannenbaum  and  collaborators  have  produced 
22  refereed  journal  papers  which  have  appeared  or  been  accepted  for  publi¬ 
cation  (papers  1-22),  4  refereed  papers  which  will  appear  as  a  book  chapters 
(papers  23-26),  14  conference  papers  (27-40).  All  of  the  conference  papers 
were  either  refereed  or  invited.  One  paper  (paper  41)  has  appeared  in  which 
Tannenbaum’s  Ph.  D.  student  supported  by  AFOSR-88-0020,  Hitay  Ozbay, 
appears  as  sole  author.  Finally  a  monograph,  Feedback  Control  Theory  (writ¬ 
ten  in  collaboration  with  John  Doyle  and  Bruce  Francis,  and  to  be  published 
by  MacMillan)  on  robust  control  was  begun  under  AFOSR-88-0020  and  will 
be  completed  under  AFOSR-90-0024. 

A  complete  list  of  these  papers  (including  the  monograph)  is  given  in 
Section  7.  This  constitutes  the  list  of  papers  produced  under  AFOSR-88- 
0020  by  Tannenbaum  and  various  collaborators.  The  grant  ended  March  31, 
1990  after  two  and  a  half  years  of  funding.  The  research  program  is  being 
continued  under  AFOSR-90-0024. 

Several  of  the  later  papers  were  also  written  during  the  intersection  of 
AFOSR-88-0020  and  the  new  contract  AFOSR-90-0024,  and  so  both  con¬ 
tract  numbers  appear  on  these  papers. 

Hitay  Ozbay  completed  his  Ph.D.  in  August  1989.  Since  then  two  grad¬ 
uate  students  Brian  Schipper  (masters)  and  Evgeny  Berzon  (Ph.D.)  have 
been  supported  on  the  contract  on  a  half  time  basis.  Brian  Schipper  will 
complete  his  masters  degree  work  in  the  spring  of  1990,  and  then  Berzon 
will  receive  full  time  support  on  the  new  contract  AFOSR-90-0024. 

We  now  sketch  the  work  done  by  A.  Tannenbaum  and  collaborators  on 
this  contract.  All  the  papers  which  have  not  been  sent  in  the  past  will  be 
sent  to  AFOSR  under  separate  cover.  (Many  of  the  papers  on  the  list  have 
already  been  sent  to  AFOSR  in  the  past.)  In  our  summary,  we  will  highlight 
a  few  of  the  key  results  and  concepts.  All  the  details  may  be  found  of  course 
in  the  papers  listed  in  Section  7. 

2  H°°  Control  of  Distributed  Systems 

During  this  contract,  the  H°°  problem  in  its  most  general  form  (the  four 
block  problem  or  standard  problem)  has  been  solved  for  arbitrary 
lumped  and  stable  distributed  systems  [8],  [11],  [14].  (As  this  contract  is 
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ending,  the  solution  is  being  extended  to  unstable  distributed  systems  with 
a  finite  number  of  unstable  poles.  This  work  will  be  described  in  our  yearly 
report  for  AFOSR-90-0024.) 

In  solving  such  problems,  we  have  introduced  a  new  class  of  opera¬ 
tors  called  skew  Toeplitz,  which  seem  ideally  suited  for  studying  H 00 
design  problems,  especially  for  distributed  parameter  systems  [1],  [3],  [5], 
[7],  [13].  These  operator  theoretic  methods  are  quite  powerful  since  they 
allow  one  to  synthesize  optimal  and  sub-optimal  compensators  just  using 
the  input/output  operators.  An  important  new  feature  of  this  approach  is 
that  the  complexity  of  the  computations  only  depends  on  the  the  weighting 
matrices  (modelling  the  disturbances)  and  not  on  the  plant  (which  may  be 
distributed).  Since  the  weighting  matrices  are  typically  taken  to  be  rational 
of  low  order,  this  approach  seems  to  fit  in  very  well  even  in  the  finite  dimen¬ 
sional  case  for  plants  with  large  state  spaces  (e.g.,  finite  element  models  of 
large  space  flexible  structures.).  Computer  implementations  for  our  proce¬ 
dures  have  been  worked  out  at  at  the  Systems  Research  Center  of  Honeywell, 
Minneapolis,  and  the  University  of  Minnesota.  Actually,  for  problems  even 
involving  finite  dimensional  systems  our  algorithm  seemed  to  converge  much 
more  quickly  than  competing  computational  schemes. 

We  now  would  like  to  describe  briefly  some  of  the  ideas  behind  our 
approach.  Full  details  can  of  course  be  found  in  the  references  given  above. 
We  want  to  indicate  how  operator  theoretic  methods  may  be  employed  to 
solve  a  very  general  case  of  the  standard  or  four  block  problem  in  H°° 
design  valid  for  a  large  class  of  distributed,  i.e.,  infinite  dimensional  systems. 

The  motivations  for  studying  the  H°°  optimization  in  systems  theory  lie 
in  the  most  natural  problems  of  control  engineering  such  as  robust  stabiliza¬ 
tion,  sensitivity  minimization,  and  model  matching.  It  can  be  shown  that, 
in  the  sense  of  H°°  optimality,  these  problems  are  equivalent,  and  can  be 
formulated  as  one  standard  problem. 

Now  it  is  quite  well-known  that  an  optimal  solution  of  the  standard 
problem  can  be  reduced  to  finding  the  singular  values  and  vectors  of  a  cer¬ 
tain  operator  (the  so-called  four  block  operator)  which  will  be  defined 
below.  Depending  on  the  specific  problem  considered,  the  corresponding 
four  block  operator  can  be  simplified  to  a  2-block  (mixed  sensitivity)  or  a 
1-block  (sensitivity)  operator.  * 

Besides  appearing  in  the  most  general  H°°  synthesis  problems,  the  four 
block  operators  also  have  a  number  of  intriguing  mathematical  properties  in 
the  sense  that  they  are  natural  extensions  of  both  the  Hankel  and  Toeplitz 
operators.  For  this  reason  they  fit  into  the  skew  Toeplitz  framework  devel- 
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oped  in  [5].  Once  again  for  the  full  details  of  our  arguments  and  details 
about  the  skew  Toeplitz  theory  applied  to  this  problem  we  refer  the  reader 
to  above  references.  Here  we  will  just  consider  the  four  block  problem  for 
single  input  /  single  output  systems. 

More  precisely,  invoking  the  Youla  parametrization  and  employing  stan¬ 
dard  manipulations  involving  inner-outer  factorizations,  for  a  large  class  of 
distributed  systems  we  may  reduce  the  standard  problem  mentioned  above 
to  the  following  mathematical  one.  Let  w,f,g,h,e  H°°,  where  w,f,g,h 
are  rational  and  m  is  nonconstant  inner.  (All  of  our  Hardy  spaces  will  be 
defined  on  the  unit  disc  D  in  the  standard  way.)  Set 


H  :=  inf {(( 


w  —  mq 
9 


S 

h 


[joo  :  9  €  H°°}. 


(1) 


Then  we  want  to  give  an  algorithm  for  calculating  the  quantity  n,  and  for 
finding  the  corresponding  qopi  G  H°°,  i.e.,  qopt  is  such  that 


H  := 


w  -  mqopt 
9 


f 

h 


Hoc 


Note  that  for  /  =  g  =  h  =  0,  this  reduces  to  the  classical  Nehari  problem. 

Our  program  will  be  to  identify  p  as  the  norm  of  a  certain  “four  block 
operator,”  and  then  give  a  determinants!  formula  for  its  computation. 


2.1  The  Four  Block  Operator 

We  will  now  define  the  four  block  operator.  We  set  H{m )  :=  H 2  0  mH 2 , 
L(m)  :=  L2  0  mH2,  and  we  let  PH( m)  :  H 2  — ►  H(m),  PL(m)  '■  L2  L(m) 
denote  the  corresponding  orthogonal  projections.  Let  5  :  H2  —*  H2  denote 
unilateral  shift,  T  :  H(m)  —  H{m)  the  compression  of  S,  and  let  U  :  L2  — 
L2  denote  bilateral  shift,  with  T(m) :  L(m)  —  L(m)  the  compression  of  U. 
Then  for  w,f,g,h  €  H°°  rational,  we  set 


A  := 


PHm)V(S)  Pnm)f(U) 

9(S)  h{U) 


Note  that 


w(T)PH{m)  f{T{m))PL(m) 

9(S)  KU) 
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(Clearly  A  :  H2  ©  L2  -  L(m)  ©  Z,2.) 

Using  the  commutant  lifting  theorem,  one  may  prove  the  following: 
Proposition  1  Notation  as  above.  Then  ||A||  =  /z. 

Thus  in  order  to  solve  the  four  block  problem  we  are  required  to  compute 
the  norm  of  the  operator  A. 

In  order  to  do  this,  we  will  first  need  to  identify  the  essential  norm  of 
A  (denoted  by  M||e).  cre  will  denote  the  essential  spectrum,  and  A(D)  will 
stand  for  the  set  of  analytic  functions  on  D  which  are  continuous  on  the 
closed  disc  D.  We  can  now  state  the  following  result  [11]: 

Theorem  1  Notation  as  above.  Let  w,f,g,h  6  A(D),  and  set 


a  max{|| 

w(C) 

.  9(C) 

/(C) 

h(C) 

II  :  C  6  °e(T)} 

(2) 

) 3  :=  max{|| 

0 

. 9(0 

0 

HO 

1 

W-CedD} 

(3) 

7  :=  sup{|| 

/(O' 
h(  C). 

ll;C€ 

dD}. 

(4) 

Then 

MIU  =  max(o,/J,7).  (5) 

2.2  Singular  System 

In  this  section,  we  will  study  the  invertibility  of  certain  skew  Toeplitz  oper¬ 
ators  [5]  which  occur  as  the  basic  elements  in  our  procedure  for  computing 
the  norm  and  singular  values  of  the  four  block  operator.  We  will  show  that 
the  calculation  of  the  singular  values  of  the  four  block  operator  A  amounts 
to  inverting  two  ordinary  Toeplitz  operators,  and  essentially  inverting  an 
associated  skew  Toeplitz  operator.  The  Fredholm  conditions  on  the  invert¬ 
ibility  of  the  skew  Toeplitz  operator  (which  is  essentially  invertible),  and  the 
coupling  between  the  various  systems  (expressed  as  matching  conditions) 
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constitutes  a  certain  linear  system  of  equations  called  the  singular  system 
which  allows  one  to  determine  the  invertibility  of  A. 

Let  p  >  max(a,/?,7).  Note  that  when  ||,4[|  >  ||A||e,  ||A||2  is  an  eigenvalue 
of  AA*.  By  slight  abuse  of  notation,  £  will  denote  a  complex  variable  as 
well  as  an  element  of  dD  (the  unit  circle).  The  context  will  always  make 
the  meaning  clear.  Of  course,  if  (  E  dD,  then  (  =  1/C- 

We  take  w,f,g,h  to  be  rational,  and  so  we  can  express  w  =  a/q,  f  =  b/q, 
g  =  c/q,  h  =  d/q ,  where  a,b,c,d,q  are  polynomials  of  degree  <  n.  Then  we 
have  that 

JWf)(S)  JWfXtO 


A  := 


(f)(5)  (i)(U) 


Now  p2  is  an  eigenvalue  of  AA *  if  and  only  if 


'  p2q(T(m))q(T(m)r 

0  1 

u 

0  P2q(U)q(Uy  j 

V 

'  / 

’ Hm)a(S )  PL(m)b(U)' 

a{S)mP 

pc(uy ' 

u 

c(S)  d(U) 

b(uypL(m ) 

d(uy 

V 

(6) 


=  0, 


for  some  non- zero 


u 

v 


where  P  :  L2 
Set 

and 


G  i(m)  ®  L2 
H 2  denotes  orthogonal  projection. 
u+  :=  Pu ,  u_  :=  (/  -  P)u 


v+  :=  Pv,  v _  :=  (I  -  P)v,  :=  (/  -  PH(m))v. 

Then  we  can  write  (6)  equivalently  as 

p2q(T(m))q(T(m)r  -  6(T(m))6(r(m))*  -b(T(m))PL(m)d(Ur 
-d(U)b(T(m)Y 


a{T)a(T)m  a(T)PH(m)c(S)' 
c(S)a(T)‘  c(S)c(ST 


p2q(U)q(Ur  -  d(U)d(Uy 


=  0 


u+ 

Let  V  :=  U*\L2  0  H2 .  If  we  apply  (I  -  P)  to  both  rows  of  (7),  we  see  that 
the  basic  block  operator  applied  to 


u_ 

u_ 


u 

V 


(7) 
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is 


'  p2q(V)q.{V)  -  b(V)b.(V)  -b(V)d.(V) 

-d(V’)b.(V)  P*q(V')q.(V)  -  d{V)d.{V) 

u_ 

v_ 

Next  applying  (I-Pn(m))  to  both  rows  of  (7),  we  see  that  the  basic  operator 
applied  to  v++  is 

C++7m>++  :=  P{(p2\q\2  -  \c\2  -  |d|2)mi;++}.  (9) 


Finally,  applying  Pn(m)  to  (7),  we  derive  that  the  basic  operator  applied  to 


u+ 

PH[m)V+ 


is 


C+ 


«+ 

PH(m)v+ 


(10) 


p2q(T)q(TY  -  b(T)b(T)m  -  a(T)a(T)'  -b(T)PH(m)d(S)' 

-d{T)b(TY  P2q{T)q{TY  -  d(T)d(TY  -  c(T)c(T)' 


u+ 

PH(m)v+ 

The  operators  C_,C++,C+  are  ail  skew  Toeplitz  (see  [5]  for  the  precise 
definition).  In  [7],  [8],  [11],  [14]  we  show  how  to  invert  C_  and  C++  under 
the  assumption  p  >  N4||e,  and  how  to  essentially  invert  C+. 

Based  on  the  inversion  formulae  of  C_,C++,  and  the  essential  inversion 
of  the  skew  Toeplitz  operator  C+,  the  singular  system  can  be  constructed 
and  the  following  theorem  can  be  proved: 
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Theorem  2  There  exists  an  explicitly  computable  5n  X  5 n  Hermitian  ma¬ 
trix  M(p)  such  that  J5  >  max{a,/3, 7}  is  a  singular  value  of  the  four  block 
operator  A  if  and  only  if 

det  M (p)  =  0 


Remarks.  See  [11]  for  the  precise  formula  for  M(p),  and  [14],  [38]  for  the 
coresponding  formula  in  the  MIMO  case.  Notice  that  the  size  of  matrix 
M[p)  only  depends  on  the  MacMillan  degree  of  the  weighting  filters  and  not 
on  the  plant.  Hence  we  can  solve  the  H°°  problem  for  distributed  systems. 
Computer  code  for  carrying  out  the  whole  procedure  has  been  written  at 
the  Systems  Research  Center  of  Honeywell. 


2.3  Optimal  Compensators 

The  above  method  also  gives  a  way  of  computing  the  optimal  compensator 
in  a  given  four  block  problem.  Indeed,  from  the  above  determinantal  formula 
one  can  compute  the  Schmidt  pair  ip,T)  corresponding  to  the  singular  value 
s  :=  ]|A||  when  s  >  H^U,..  We  will  indicate  how  one  derives  the  optimal 
interpolant  (and  thus  the  opimal  compensator)  from  these  Schmidt  vectors. 
In  order  to  do  this,  notice 

Alp  =  ST) 

Thus,  there  exists  qopt  e  H°°  with 


(W  -  QoptWi  +  /02  =  ST)  1 


gipi  +  hip2  =  skr)2 


where 


Tp  = 


01 

77  — 

m 

02 

)  V  — 

.  % 

One  can  show,  that  ipi  ^  0,  so  that 


Qopt  —  w  — 


sm  -  fip? 

*P\ 


Note  from  qopt,  using  the  Youla  parametrization,  we  can  derive  the  cor¬ 
responding  optimal  controller  in  a  given  systems  design  problem.  See  also 
[25]  for  an  extension  of  the  theory  of  Adamjan-Arov-Krein  (valid  for  the 
Hankel  operator)  to  the  singular  values  of  the  four  block  operator  and  their 
relationship  to  more  general  interpolation  and  distance  problems. 
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2.4  Suboptimal  Compensators 

We  would  like  to  describe  now  our  techniques  for  constructing  suboptimal 
compensators  for  distributed  plants.  Indeed,  using  the  one-step  dilation 
procedure,  we  have  given  a  way  of  explicitly  parametrizing  all  the  suboptimal 
controllers  for  such  generalized  interpolation  problems  [10].  We  have  also 
studied  in  this  context  finite  dimensional  suboptimal  controllers  for  infinite 
dimensional  systems  in  [18],  [34],  [39].  Such  techniques  are  of  course  very 
important  in  the  practical  implementation  of  distributed  controllers. 

For  simplicity,  we  will  consider,  the  sensitivity  minimization  problem. 
Thus  we  are  required  to  find  to  find  an  internally  stabilizing  controller  C 
such  that  the  following  optimum  performance  is  achieved 


inf  \\W(1  +  PC) 

C  stabilizing 


-ll 


=  :  ft. 


where  P  is  the  plant  to  be  controlled  and  W  is  the  weight  modelling  the 
disturbances.  Assuming  that  the  weight  is  an  outer  function  and  the  plant 
P  is  stable  we  can  transform  this  problem  to  a  Nehari  problem  (in  the 
usual  way  of  first  invoking  the  Youla  parametrization  for  the  controller  C  = 
Qc(l  -  PQC)~X,  Qc  G  H°°,  (1  -  PQC )  ^  0;  and  then  finding  an  inner/outer 
factorization  for  P  =  mP0,  where  m  is  inner  and  PQ  is  outer): 


P  =  WW  ~  mGIU- 

Now  given  a  tolerance  e  >  0,  we  say  that  Ct  is  suboptimal ,  (or  approx¬ 
imately  optimal ),  with  tolerance  e,  if  it  internally  stabilizes  the  system  and 
satisfies  the  bound 


||W(1  +  PC()~X\\oo  <  p  +  (  =■  P- 

We  are  thus  led  to  consider  the  following  problem:  given  p  >  p,  find  the  sei 
of  all  Q  6  H 00  such  that 

\\W  -  mQHoo  <  p.  (11) 

(Notice  that  from  such  a  suboptimal  interpolant  Q,  we  can  solve  for  the 
corresponding  suboptimal  controller.) 

In  [10],  we  give  a  parametrization  of  all  such  suboptimal  interpolants. 
Indeed,  suppose  that  the  weight  is  rational:  W(z)  =  p(z)/q(z)  where  p(z)  = 
po  +  zpi  +  ...  +  znpn  and  q(z)  =  qQ  +  zqi  +  . . .  +  znqn,  (i.e.,  n  is  the 
maximum  of  the  degrees  of  p  and  q,  so  some  of  the  above  coefficients  may 
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well  be  zero).  Let  S  denote  the  unilateral  shift  on  H 2  and  define  the  space 
//(to)  =  H2  0  mH2.  Then  the  compressed  shift  associated  with  //(to)  is 
defined  as  T  :=  PH(m)S\H(m),  where  Pn{m)  denotes  orthogonal  projection. 

First,  consider  the  optimal  case:  p  =  p.  The  optimal  interpolant  Qopt, 
which  makes  |[Z?opt||oc  =  P<  where 

Popt  =  11  —  mQopt, 

can  be  computed  using  Sarason’s  theorem  which  states  that 
p  =  \\W(T)l  W{T)  :=  p(T)q(T)~l . 

The  essential  norm  can  be  defined  as  follows: 


||lF(T)||e  =  sup{|lF(£)|  :  £  singular  point  of  m}. 

Assuming  that  p  >  ||lT(T)||e,  we  have  that  \V(T)  attains  its  norm  at  a 
singular  value  p  =  p.  In  this  case  there  exists  a  singular  vector  h0,  for  the 
skew  Toeplitz  operator 

Ap  := p2q(T)q(Tr~p(T)p(T)’ 

(*  denotes  adjoint)  which  makes 


Aph0  =  0. 


The  vector  hQ  can  be  computed  explicitly  from  the  problem  data  W  =  p/q 
and  m  in  terms  of  a  determinantal  formula;  see  [7],  [10],  and  [18].  Then  Bopt 
can  be  found  via  Sarason’s  result  as 


Bopt  —  p 


q(Trho 

p(T)mh0 


Let  us  now  consider  the  case  where:  p  >  p.  It  is  obvious  that  in  this 
case  Ap  is  invertible  and  its  inverse  can  be  computed  explicitly;  again,  the 
formula  is  given  in  [10].  This  is  going  to  be  used  in  the  characterization  of 
all  the  suboptimal  solutions  Q3  £  H°°  which  make 


II W  -  mQs |U  <  P- 


(12) 


This  characterization  is  obtained  using  the  one  step  extension  procedure 
of  Adamjan- Arov-Krein.  Here  we  want  to  summarize  the  method  briefly. 
Set  mv(z)  :=  zm(z)  and  let  Tv  denote  the  compression  of  S  to  H{mv)  = 
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H(m)  ®  Cz.  For  a  £  C  fixed,  the  problem  of  finding  Bopt{z,a)  =  ( W  - 
am  -  mvQ°pt)(z)  such  that 

||B0pf(.,a)||00  =  ||(ty-om)(rt,)||  =  p 

can  be  solved  using  the  technique  described  above  for  the  optimal  case.  From 
one  step  extension  theory,  we  know  that  the  set  of  all  such  a  £  C  form  a 
circle,  say  T.  Furthermore,  the  equation  of  T  can  be  explicitly  calculated. 
Then  the  set  of  all  suboptimal  solutions  Qs  £  H°°  satisfying  (12)  is  obtained 
in  terms  of  Bopt(z,<j>(u)): 

W  -  mQa  =  Bopt(z,4>{u)), 

where  4>(  z)  is  a  linear  fractional  map  taking  the  unit  circle  to  T,  and  u  £  H°°, 
IMloo  <  1  is  the  free  parameter.  The  explicit  characterization  is  as  follows. 
Set 

9i  ■=  (. p2q(T)PH(m)q(Sr  -  p(T)PH{m]p(Sr)m , 

9i  :=  QqP(T)(1  -  m  m(0)), 


hi:=A;lgu  h2  :=  A~lg2- 


For  a  given  a  £  T  define 


ha(z)  :=  m(z)  -  hx{z)  -  &h2(z), 


B(z,a )  := 

Then  we  have  the  following  result. 


P2q(S)'hQ 
p(S)’ha  -  7Sq0' 


Theorem  3  The  set  of  all  functions  of  the  form 

B(z)  =  W{z)  -  m(z)Qs(z) 
unth  Qs  £  H°°,  such  that  U^Hoo  <  p,  is  given  by 

P\j)  ha  ago 

where  r  and  rj  are  certain  explicitly  computable  constants.  (See  [10]  and  [18] 
for  the  formulae.) 


From  the  above  parametrization  can  obtain  the  structure  of  all  subop- 
timal  H°°  controllers.  Using  the  notation  of  Theorem  3,  we  set  Ba(z)  := 
B(z,a).  We  can  find  the  controller  from  C  =  Qc(  1  -  PQc)~l ,  the  Youla 
parametrization,  where  Qc  is  such  that 


Ba  =  W  -  PWQC. 


Therefore, 

C  =  P~X(B~XW  -  1). 

We  now  study  Ba , 

B  <z\-  P2q(*)M*)~  P2K(Z) 

a  ~  p(z)ha(z)  -  hp(z )  -  q0znT$ 

where  hq(z)  and  hp(z)  are  polynomials  of  degree  <  n  —  1  and  q(z)  = 
znq(z~x),  similarly  p(z)  =  znp(z~x).  Then, 

_1_  / p(z)ha(z)  -  hp(z)  -  q0zn7S  p(z)  _  \ 
m\  p2q(z)ha(z)  -  p2hq{z)  q(z)  ) 


...  J_  ( -A(z)/t0(z)  -  p(z)hp(z)  +  p2g(z)hq(z)  -  q0znp(zp\ 

m  V  iP2q{z)ha{z)  ~  P2hq(z))q(z)  )  ' 

where  A (z)  =  p2q(z)q(z)  -  p(z)p(z).  Recall  that  ha(z )  =  m(z)  -  hi(z )  - 
7Zh,2(z).  It  is  easy  to  see  from  the  inversion  of  the  skew  Toeplitz  operator 
Ap,  that  hi  and  h2  have  the  following  form: 


and 


hi(z)  = 


h2(z)  = 


fi(z)  +  m(z)Fi(z) 
\(z) 

f2(z)  +  m(z)F2{z ) 
A(z) 


for  some  f\,Fi,f2,F2  polynomials  of  degree  <  2n.  This  leads  us  to  the 
following  expression: 

-A(z)  >>  Gu(z) 


PaC 


=  (  -*(«)  ^  _ 

\P2q(z)q(z)J  1 


where 


Gu(z)  =  q{z) 


p2q{z)q{z)J  1  +  m(z)Gu(z) 
Fa{z)  -  A(z) 


q(z)fa(.z)  +  hq(z)X(z) 
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Fa(z)  :=  Fi{z)  +  75 F2(z),  and  fa(z)  :=  /i(z)  +  75/2(z). 

and 

r 

75  =  -  +  ??. 
u 

Note  that 

-A(z)  p(z)p(z)  _  W(z)W(z -1) 
p2q(z)q(z)  p2q(z)q(z)  p2 

We  summarize  the  above  formulae  with  the  following. 


Corollary  1  The  set  of  all  controllers  which  internally  stabilize  the  plant 
P,  and  satisfy  the  bound 


\\W{\  +  PC)^U<P 


for  p  >  p,  have  the  form 

CM  = 


-  1 


Gu(z) 


1  +  m(z)Gu(z) 


Po'(z) 


«  e  H°°,  HttHoo  <  l,  where  Gu(z )  is  a  linear  fractional  transformation  in 
the  free  parameter  u: 


/-/.x  <Pi  (*)  +  ¥*(*)« 

\  /  s. 

<pi ,<f4  polynomials  of  degree  <  3  n.  They  can  be  computed  explicitly 
from  the  equations  given  in  [10]  and  [18]  via  fi,  F\,  f2,  F2,r  and  tj.  □ 

In  [18],  [34],  [39],  we  discuss  how  to  employ  these  methods  in  finding 
finite  dimensional  suboptimal  controllers  for  distributed  systems.  In  partic¬ 
ular,  in  [34]  these  techniques  are  used  for  the  mixed  sensitivity  design  of  a 
flexible  beam  modelled  by  the  Euler- Bernoulli  equation.  (This  work  is  being 
continued  under  AFOSR-90-0024,  and  a  full  journal  paper  is  being  prepared 
on  this  subject.) 
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2.5  Other  Results  in  H°°  Control 

We  have  also  verified  a  formula  which  combines  the  state  and  frequency 
approaches  to  /f°°-optimal  design  in  the  one  block  (sensitivity  minimiza¬ 
tion)  case  [2].  Indeed,  the  formula  utilizes  the  fact  that  the  weighting  filter 
is  taken  to  be  rational  while  the  plant  may  be  distributed.  Work  on  our 
new  contract  AFOSR-90-0024  is  being  continued  to  extend  this  to  the  more 
general  two  and  four  block  frameworks.  This  type  of  formula  combines  the 
advantages  of  both  the  state  space  (Ricatti  equation)  and  frequency  domain 
(input/output)  approaches  to  H°°  theory. 

Moreover  in  [25],  we  have  given  a  precise  interpolation  theoretic  descrip¬ 
tion  of  the  singular  values  of  the  four  block  operator  generalizing  the  classical 
results  of  Adamjan-Arov-Krein.  We  believe  this  type  of  result  should  prove 
very  useful  in  approximation  theory  (along  the  lines  of  the  present  applica¬ 
tions  of  the  Adamjan-Arov-Krein  result). 

A  general  description  of  our  interpolation  operator  theoretic  approach 
to  distributed  parameter  systems  is  given  in  [4]  and  [24].  A  monograph  be¬ 
ing  written  jointly  by  J.  Doyle,  B.  Francis,  and  A.  Tannenbaum  (which  was 
started  under  AFOSR-88-0020  and  which  will  be  completed  under  AFOSR- 
90-0024),  Feedback  Control  Theory  (to  be  published  by  MacMillan),  de¬ 
scribes  from  a  basic  point  of  view  some  of  the  uses  of  interpolation  techniques 
in  robust  control  theory. 

3  Nonlinear  Systems 

A  large  portion  of  the  work  on  AFOSR-88-0020  has  been  devoted  to  robust 
nonlinear  design.  This  work  is  being  carried  over  into  AFOSR-90-0024  as 
well.  We  believe  that  we  now  have  a  viable  nonlinear  generalization  of  the 
(linear)  H°°  theory  in  both  the  weighted  sensitivity  minimization  (one  block) 
and  mixed  sensitivity  (two  block)  frameworks.  This  is  based  on  a  iterative 
commutant  lifting  theorem  [12],  [31],  [35]  which  gives  an  explicit  design 
procedure  for  nonlinear  systems  and  captures  the  tf^-control  problem  for  a 
large  class  of  nonlinear  plants.  We  have  also  defined  a  notion  of  rationality 
for  nonlinear  systems,  and  we  have  proven  that  the  iterative  commutant 
lifting  procedure  produces  rational  controllers  (in  this  nonlinear  sense)  if  we 
start  from  rational  data  [9],  We  have  thus  been  able  to  write  computer  code 
for  this  procedure  along  the  lines  that  was  done  for  the  four  block  problem 
using  the  theory  of  skew  Toeplitz  operators. 

We  would  like  to  explain  a  bit  our  approach  to  nonlinear  H°° .  We  will 
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concentrate  below  on  the  one  block  problem.  Extensions  have  already  gone 
through  to  the  mixed  sensitivity  (two  block)  problem. 

3.1  Analytic  Input/Output  Mappings 

In  order  to  carry  out  our  extension  of  H°°  synthesis  theory  to  nonlinear 
systems,  we  will  need  to  first  discuss  a  few  standard  results  about  analytic 
mappings  on  Hilbert  spaces. 

Let  G  and  H  denote  complex  Hilbert  spaces.  Set 

Bro(G)  :={geG:  |!s||  <  r0} 

(the  open  ball  of  radius  r0  in  G  about  the  origin).  Then  we  say  that  a 
mapping  <p :  Bro(G )  — ►  H  is  analytic  if  the  complex  function  (zi, . . . ,  zn)  t-> 
((•*i5i  +  . . .  +  zngn),h)  is  analytic  in  a  neighborhood  of  (1, 1, . . . ,  1)  E  Cn 
as  a  function  of  the  complex  variables  zi,...  ,zn  for  all  glt . . .  ,gn  6  G  such 
that  ||<7!  +  . . .  +  <7n||  <  r0,  for  all  h  €  H,  and  for  all  n  >  0.  (Note  that  we 
denote  the  Hilbert  space  norms  in  G  and  H  by  ||  ||  and  the  inner  products 

by  {,)•) 

We  will  now  assume  that  <p(0)  =  0.  It  is  easy  to  see  that  if  cp  :  BT„(G)  — ► 
H  is  analytic,  then  <p  admits  a  convergent  Taylor  series  expans  ion,  i.e. 

4>{g)  =  Ms)  +  Ms, s)  +  •  •  •  +  Ms,  ••*,?)  +  •■• 

where  <pn  :  G  x  •  •  •  x  G  — ►  H  is  an  n-linear  map. 

Clearly,  without  loss  of  generality  we  may  assume  that  the  n-linear  map 
ls\,  •  •  •  ,9n)  Ml,-  •  •  ,9n)  is  symmetric  in  the  arguments  gi,---,9n-  For 
<t>  a  Volterra  series,  <pn  is  basically  the  n</l-Volterra  kernel. 

Now  set 

Mgi  ®---®gn)  ••=  Msu---,Sn)- 

Then  <j>n  extends  in  a  unique  manner  to  a  dense  set  of  G®n  :=  G  ®  ®  G 

(tensor  product  taken  n  times).  Notice  by  G®n  we  mean  the  Hilbert  space 
completion  of  the  algebraic  tensor  product  of  the  G’s.  Clearly  if  <pn  has 
finite  norm  on  this  dense  set,  then  4>n  extends  by  continuity  to  a  bounded 
linear  operator  <pn  :  G®n  — ►  H.  By  abuse  of  notation,  we  will  set  <pn  :=  <pn. 
We  now  make  the  following  definitions. 

Definitions  1. 

(i)  Notation  as  above.  By  a  majorizing  sequence  for  the  holomorphic  map  d>, 
we  mean  a  positive  sequence  of  numbers  a  „  n  =  1,2,.. .  such  that  ||d>n||  <  an 
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for  n  >  1.  Suppose  that  p  :=  limsupan,1/**  <  oo.  Then  it  is  completely 
standard  that  the  Taylor  series  expansion  of  <f>  converges  at  least  on  the  ball 
Br(G)  of  radius  r  =  1/p. 

(ii)  If  4>  admits  a  majorizing  sequence  as  in  (i),  then  we  will  say  that  (p  is 
majorizable. 

3.2  Nonlinear  Control  Problem 

We  will  describe  in  this  section  the  control  problem  in  which  we  are  in¬ 
terested.  First,  we  will  need  to  consider  the  precise  kind  of  input/output 
operator  we  will  be  considering. 

We  assume  that  all  of  the  operators  we  consider  are  causal  and  majoriz¬ 
able.  H2(Ck)  denotes  the  standard  Hardy  space  of  Ck-valued  functions  on 
the  unit  circle  ( k  may  be  infinite,  i.e.,  in  this  case  Ck  is  replaced  by  h2,  the 
space  of  one-sided  square  summable  sequences).  We  now  make  the  following 
definition: 

Definition  2. 

Let  5  :  H2( Ck)  — ►  H2(Ck)  denote  the  canonical  unilateral  right  shift. 
Then  we  say  an  input/output  operator  <p  is  locally  stable  if  it  is  causal 
and  majorizable,  <£(0)  =  0,  and  if  there  exists  an  r  >  0  such  that  <p  : 
Bt{H2{ Ck))  -  H2{ Ck)  with  S<f>  =  <f>o  S  on  Br(H2(Ck)).  We  set 

Ci  :=  {space  of  locally  stable  operators}. 

Since  the  theory  we  are  considering  is  local,  the  notion  of  local  stability  is 
sufficient  for  all  of  the  applications  we  have  in  mind. 

The  theory  we  are  about  to  give  holds  for  all  plants  which  admit  coprime 
locally  stable  factorizations.  However,  for  simplicity  we  will  assume  that 
our  plant  is  also  locally  stable.  Accordingly,  let  P,  W  denote  locally  stable 
operators,  with  W  invertible.  As  before,  P  represents  the  plant,  and  W  the 
weight  or  filter  on  the  set  of  disturbances  whose  energy  is  bounded  by  a  fixed 
constant.  Now  we  say  that  the  feedback  compensator  C  locally  stabilizes  the 
closed  loop  if  the  operators  (I  -f-  P  o  C)~l  and  C  o  (/  -(-  P  o  C)-1  are  well- 
defined  and  locally  stable.  One  can  show  that  C  locally  stabilizes  the  closed 
loop  if  and  only  if 

C  =  qo(I-  Poq)~l  (13) 
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for  some  q  €  C/.  Notice  then  that  the  weighted  sensitivity  (/  +  P  o  C)~l  o  W 
can  be  written  as  W-Poq,  where  q  :=  qoW .  (Since  W  is  invertible,  the  data 
q  and  q  are  equivalent.)  In  this  context,  we  will  call  such  a  q,  a  compensating 
parameter.  Note  that  from  the  compensating  parameter  q,  we  get  a  locally 
stabilizing  compensator  C  via  the  formula  (13). 

The  problem  we  would  like  to  solve  here,  is  a  version  of  the  classical 
disturbance  attenuation  problem.  This  of  course  corresponds  to  the  “mini¬ 
mization”  of  the  “sensitivity”  W  -  P  o  q  taken  over  all  locally  stable  q.  In 
order  to  formulate  a  precise  mathematical  problem,  we  need  to  say  in  what 
sense  we  want  to  minimize  W-Poq.  This  we  will  do  in  the  next  section 
where  we  will  propose  a  notion  of  “sensitivity  minimization”  which  we  seems 
quite  natural  to  analytic  input/output  operators. 

3.3  Nonlinear  Sensitivity  Function 

We  now  define  a  fundamental  object,  namely  a  nonlinear  version  of  sensi¬ 
tivity.  We  will  see  that  while  the  optimal  H°°  sensitivity  is  a  real  number  in 
the  linear  case,  the  measure  of  performance  which  seems  to  be  more  natural 
in  this  nonlinear  setting  is  a  certain  function  defined  in  a  real  interval.  This 
new  kind  of  performance  criterion  is  one  of  the  keys  concepts  developed  in 
[9]  and  [12]. 

In  order  to  define  our  notion  of  sensitivity,  we  will  first  have  to  partially 
order  germs  of  analytic  mappings.  All  of  the  input/output  operators  here 
will  be  locally  stable.  We  also  follow  here  our  convention  that  for  given  4>  £ 
Ci,  <t>n  will  denote  the  bounded  linear  map  on  the  tensor  space  {H2( Ck))®n 
associated  to  the  n-linear  part  of  which  we  also  denote  by  </>n  (and  which 
we  always  assume  without  loss  of  generality  is  symmetric  in  its  arguments). 
The  context  will  always  make  the  meaning  of  <pn  clear. 

We  can  now  state  the  following  definitions: 

Definitions  3. 

(i)  For  W,P,q  £  C\  (W  is  the  weight,  P  the  plant,  and  q  the  compensating 
parameter),  we  define  the  sensitivity  function  S(q), 

S(q)(p)  :=  XyiKW-Po^H 

n=l 

for  all  p  >  0  such  that  the  sum  converges.  Notice  that  for  fixed  P  and  W, 
for  each  q  £  Ci,  we  get  an  associated  sensitivity  function. 
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(ii)  We  write  S(q)  ■<  S(q ),  if  there  exists  a  p0  >  0  such  that  S(q)(p)  < 
S(q)(p)  for  all  p  g  [0,po].  If  S(q )  <  S(q)  and  S(q)  <  S(q),  we  write 
S(q)  S  S(q).  This  means  that  S{q){p)  =  S(q)(p)  for  all  p  >  0  sufficiently 
small,  i.e.,  S(q)  and  S(q)  are  equal  as  germs  of  functions. 

(iii)  If  S(q)  <  S(q),  but  S(q )  2?  £(9),  we  will  say  that  q  ameliorates  9.  Note 
that  this  means  S(q)(p)  <  S(q)(p)  for  all  p  >  0  sufficiently  small. 

Now  with  Definitions  3,  we  can  define  a  notion  of  “optimality”  relative 
to  the  sensitivity  function: 

Definitions  4. 

(i)  q0  €  Ci  is  called  optimal  if  S(q0)  <  S(q)  for  all  q  6  C/. 

(ii)  We  say  q  €  Ci  is  optimal  with  respect  to  its  n-th  term  qn,  if  for  every 
n-linear  qn  €  Ci ,  we  have 

+  •  •  ■  +  9n—  1  +  9n  +  9n+ 1  •  ••)  ^  +  •  •  •  +  9n- 1  +  9n  +  9n+ 1  +  •••)• 

If  q  6  Ci  is  optimal  with  respect  to  all  of  its  terms,  then  we  say  that  it  is 
partially  optimal. 

3.4  Iterative  Commutant  Lifting  Method 

In  this  section,  we  discuss  the  main  construction  of  this  paper  from  which  we 
will  derive  both  partially  optimal  and  optimal  compensators  relative  to  the 
sensitivity  function  given  in  Definitions  3  above.  As  before,  P  will  denote  the 
plant,  and  W  the  weighting  operator,  both  of  which  we  assume  are  locally 
stable.  As  in  the  linear  case,  we  always  suppose  that  Pi  is  an  isometry,  i.e. 
Pi  is  inner. 

We  begin  by  noting  the  following  key  relationship: 

( W-Poq)k  =  Wk -  ^  E 

l<j<k  ii  +....+«y=fc 

Note  that  once  again  for  <j>  locally  stable,  </>„  denotes  the  n-linear  part  of  <fi, 
as  well  as  the  associated  linear  operator  on  the  appropriate  tensor  space. 

We  are  now  ready  to  formulate  the  iterative  commutant  lifting  procedure. 
Let  II  :  H2( Ck)  — +  H2{ Ck)  Q  Pi/f2(Ck)  denote  orthogonal  projection. 
Using  the  linear  commutant  lifting  theorem  (CLT),  we  may  choose  qi  such 
that 

||Wi  -  Pifill  =  linwi||. 
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Now  given  this  qly  we  choose  (using  CLT)  q2  such  that 

\\w2  -  p2(qi  ®  qi)  -  piq2 1|  =  ||n (w2  -  p2(qi  ®  9l))||. 

Inductively,  given  q\, . .  .,qn- 1,  set 

An  :=  (W„  -  H  PMh  )) 

for  n  >  2.  Then  from  the  CLT,  we  may  choose  qn  such  that 

\\An-piqn\\  =  ||im„||.  (14) 

We  now  come  to  the  key  point  on  the  convergence  of  the  iterative  commutant 
lifting  method. 

Proposition  2  With  the  above  notation,  let  q'A)  :=  qi  +  q2  + _  Then 

qW  e  C,. 

Note  that  given  any  q  6  C/,  we  can  apply  the  iterative  commutant  lifting 
procedure  to  W  -  P  o  q.  Now  set 

5n(q)(p):=E^lln(^-Po^H- 

n=l 

Clearly,  Sn(<?)  <  S(q)  (as  functions).  We  can  now  state  the  following  result 
whose  proof  is  immediate  from  the  above  discussion: 

Proposition  3  Given  q  6  Ci,  there  exists  q  6  Cj,  such  that  S(q)  =  5n(g). 
Moreover  q  may  be  constructed  from  the  iterated  commutant  lifting  proce¬ 
dure. 

Moreover,  we  easily  have  the  following  result: 

Proposition  4  q  is  partially  optimal  if  and  only  if  S(q)  ^  Sn(<?). 

Finally,  from  the  above  discussion,  we  can  prove: 
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Theorem  4  For  given  P  and  W  as  above,  any  q  £  C/  is  either  partially 
optimal  or  can  be  ameliorated  by  a  partially  optimal  compensating  parameter. 

It  is  important  to  emphasize  that  a  partially  optimal  compensating  pa¬ 
rameter  need  not  be  optimal  in  the  sense  of  Definition  4(i).  Basically,  what 
we  have  shown  here  is  that  using  the  iterated  commutant  lifting  procedure, 
we  can  ameliorate  any  given  design.  The  question  of  optimality  will  be 
considered  next. 

3.5  Optimal  Nonlinear  Compensators 

We  will  mention  now  some  of  our  results  about  optimal  compensators.  Ba¬ 
sically,  we  want  to  give  conditions  when  the  the  iterated  commutant  lifting 
procedure  leads  to  an  optimal  design.  (We  have  seen  that  it  always  leads 
to  a  partially  optimal  one.)  We  begin  with  the  following  very  general  result 
from  [12]: 


Theorem  5  There  always  exist  optimal  compensators. 


In  fact,  in  several  very  important  cases,  we  can  use  the  iterated  com¬ 
mutant  lifting  procedure  to  explicitly  construct  optimal  compensators.  We 
now  quote  one  of  these  cases  from  [12]: 


Theorem  6  Let  P  and  W  be  SISO,  locally  stable,  with  the  linear  part  of  P 
rational.  Then  the  partially  optimal  compensating  parameter  qopt  constructed 
by  the  iterated  commutant  lifting  procedure  is  optimal. 


In  collaboration  with  colleagues  at  Honeywell,  we  have  been  applying  these 
methods  to  certain  systems  with  input  saturations  (e.g.,  a  saturated  double 
integrator)  in  order  to  attenuate  bounded  energy  output  disturbances  [35]. 
Much  more  work  is  being  performed  in  this  direction  which  will  be  part  of 
AFOSR-90-0024. 
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3.6  Nonlinear  Beurling-Lax  Theory 

There  were  several  theoretical  results  in  the  area  of  nonlinear  systems  during 
the  AFOSR-88-0020  funding  which  we  mention  briefly.  These  were  based  on 
a  nonlinear  Beurling-Lax-Halmos  theorem  [6].  The  linear  Beurling- 
Lax- Halmos  theorem  characterizes  the  shift-invariant  closed  subspaces  of  a 
complex  separable  Hilbert  space  H . 

In  [6],  we  formulate  a  theory  relative  to  locally  shift-invariant  nonlinear 
analytic  manifolds  in  Hilbert  space,  and  develop  the  corresponding  repre¬ 
sentation  theory.  We  should  note  the  proof  of  the  nonlinear  Beurling-Lax- 
Halmos  theorem  makes  use  of  Poincare-Dulac  type  changes  of  variable  from 
the  theory  of  nonlinear  ordinary  differential  equations  in  the  parameter  space 
of  the  manifold  as  well  as  the  commutant  lifting  theorem. 

These  ideas  allow  us  to  give  nonlinear  analogues  of  the  classical  notion 
of  inner-outer  factorization,  and  should  be  applicable  to  certain  problems  in 
nonlinear  realization  theory,  and  a  number  of  nonlinear  optimization  prob¬ 
lems.  They  also  lead  to  a  new  linearization  result  for  very  general  nonlinear 
input /output  operators,  as  well  as  nonlinear  analogues  of  the  classical  Han- 
kel  and  Toeplitz  operators.  (In  AFOSR-90-0024,  we  will  be  exploring  the 
possibilities  of  globalizing  this  theory.) 

4  Multivariable  Gain  Margin  Problem 

We  have  found  an  analytic  procedure  for  solving  the  multivariable  gain  mar¬ 
gin  problem,  based  on  a  novel  interpolation  scheme  which  we  call  spectral 
Nevanlinna-Pick  interpolation  [19],  [21],  [22],  [27],  [37].  This  involves  ma¬ 
trix  interpolation  in  which  one  bounds  the  spectral  radius,  and  not  the 
norm  of  the  interpolants  as  in  ordinary  Nevanlinna-Pick  theory.  (Ordinary 
Nevanlinna-Pick  is  precisely  the  type  of  mathematical  problem  that  arises 
in  standard  H°°  synthesis.)  The  spectral  problem  has  been  solved  in  both 
its  matricial  and  tangential  (directional)  forms.  We  would  like  to  give  a  few 
details  about  these  ideas  now. 

A  key  observation  is  that  the  standard  one  block  H°°  optimal  control 
problem  in  the  finite  dimensional  case  amounts  to  Nevanlinna-Pick  interpo¬ 
lation  (both  for  SISO  and  MIMO  systems).  Indeed,  the  connection  between 
interpolation  theory  and  certain  questions  involving  LTI,  finite  dimensional 
plants  is  the  simple  fact  that  the  problem  of  internal  stabilization  reduces  to 
Lagrange-type  interpolation  for  such  plant.  The  basis  of  H°°  optimization 
theory  is  that  for  certain  design  problems,  one  is  required  to  interpolate 
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on  the  unit  disc  (or  equivalently,  the  right  half  plane)  by  analytic  matrix- 
valued  functions  of  bounded  norm,  which  is  precisely  the  Nevanlinna-Pick 
interpolation  problem. 

For  SISO  systems  gain  margin  optimization  works  in  the  same  way  in 
that  via  a  conformal  equivalence,  one  can  transform  the  question  of  finding 
an  interpolating  function  whose  range  is  a  certain  simply-connected  subset 
of  the  complex  plane  to  one  whose  range  is  the  unit  disc,  and  hence  once 
more  derive  an  interpolation  problem  of  the  Nevanlinna-Pick  kind.  Now  for 
MIMO  systems,  we  have  shown  that  if  one  wants  to  play  the  same  game 
with  multivariable  generalizations  of  gain  margin,  one  derives  interpolation 
not  with  a  norm  constraint  but  with  a  spectral  radius  constraint. 

4.1  Internal  Stability  and  Interpolation 

We  would  first  like  to  briefly  discuss  how  the  problem  of  internal  stabilization 
for  LTI,  finite  dimensional  plants  reduces  to  one  of  interpolation. 

Let  P(s )  denote  a  pxm  LTI  finite  dimensional  plant  and  C(s )  an  mxp 
internally  stabilizing  compensator.  In  the  usual  way  we  define  the  sensitiv¬ 
ity  function  to  be  S(s)  :=  (/  +  P(s)C(s))~l .  Then  invoking  the  standard 
coprime  factorizations  we  get  that 

S(s)  -  L  -  U Z 

where  L  and  U  are  completely  determined  by  P(s),  U  is  an  inner  matrix- 
valued  function,  and  the  “free  parameter”  Z  G  RH^Xp  (—  the  space  of  mxp 
matrices  with  entries  which  are  reol  r'tion'l  functor.:  hounded  in  the  right 
half  plane  H).  Now  one  wants  to  compute 

inf{||S|U  :  C  stabilizing  compensator}. 

As  is  well-known  this  problem  can  either  be  reduced  to  one  of  tangential 
Nevanlinna-Pick  interpolation,  or  matricial  Nevanlinna-Pick  interpolation 
via  the  tranformation 

Ual9S  =  Ual9L-{dztU)IZ, 

where  Ual 9  denotes  the  algebraic  adjoint  of  U,  and  I  the  identity  matrix. 

With  these  preliminary  remarks  made  we  are  ready  to  formulate  our  con¬ 
trol  problem  which  will  lead  to  the  spectral  Nevanlinna-Pick  interpolation. 
Let  P(s)  be  a  pxm  MIMO  finite  dimensional  plant.  Consider  the  following 
family  of  plants 

Pk(s)  :=  {kP(s)  :  k  G  K} 
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K:=  {*eC:*=l+J,M<r}. 

Then  we  may  show  that  C(s)  internally  stabilizes  the  closed  loop  for 
the  family  Pk(s),  k  6  A'  if  and  only  if  there  exists  a  rational  matrix-valued 
function  S(s)  which  is  analytic  and  bounded  in  the  right  half  plane  such 
that  5  :  H  —*  G  and  which  moreover  satisfies  standard  Nevanlinna-Pick 
interpolation  conditions,  where  H  denotes  the  closed  right  half  plane  U  {oo}, 
and 

G  {p  x  p  matrices  M  :  det(/  +  kM)  ^  0,  k  6  K}. 

But  now  it  is  easy  to  construct  a  conformal  equivalence  (which  is  a 
linear  fractional  transformation)  ip  :  G  Cl  where  Q  denotes  the  space  of 
p  x  p  matrices  with  spectral  radius  less  than  one.  From  standard  conformal 
mapping  theory,  the  interpolation  constraints  on  S  may  be  tranformed  to 
similar  constraints  on  <f>o S  :  H  — ►  ft.  In  other  word?,  we  have  a  Nevannlinna- 
Pick  type  problem  in  which  instead  of  bounding  the  norm  of  the  interpolant 
we  bound  its  spectral  radius. 

We  should  also  add  that  a  similar  trick  works  for  real  parameter  varia¬ 
tions  ( k  6  [a,h]  a  real  interval  with  0  <  a  <  1  <  6),  and  we  can  even  take 
parameter  variation  spaces  of  the  form 

I<  :=  {diag(fci,...,fcp)} 

where  the  kj's  may  be  either  real  or  complex  (p  in  the  lxl  block  case). 

4.2  Solution  of  the  Spectral  Nevanlinna-Pick  Problem 

We  now  outline  our  solution  to  the  spectral  Nevanlinna-Pick  interpolation 
problem.  For  simplicity,  we  will  consider  just  the  matricial  case.  (The  tan¬ 
gential  case  is  worked  out  in  [21].)  Following  standard  mathematical  practice 
we  work  on  the  unit  disc  D  rather  than  (its  conformal  equivalent)  the  right 
half  plane  H.  Accordingly,  we  are  given  n  distinct  points  zi,...,z„  6  D, 
and  n,  p  x  p  matrices  Pi,.. .  ,F„.  We  want  necesssary  and  sufficient  condi¬ 
tions  for  the  existence  of  a  rational  matrix-valued  function  F(z),  analytic  in 
the  unit  disc  such  that 

F(zt)  =  Ft 

for  i  ss  1, . . .  ,n,  and  such  that 

sup{sr(P(z))  :  z  €  D}  <  1,  (15) 
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where  for  a  matrix  Q,  sr(Q)  denotes  its  spectral  radius.  We  can  now  state 
the  following  theorem: 

Theorem  7  With  the  above  notation,  F(z)  exists  if  and  only  if  there  exist 
invertible  p  X  p  matrices  Mi  for  i  —  1 , . . . ,  n  such  that 

I  -  MiFiMrx{MjFjM-xy 

1  -  ZiZj 

We  would  now  like  to  discuss  a  bit  the  derivation  of  Theorem  7  which 
will  also  give  us  a  numerical  scheme  for  the  theorem’s  implementation.  Let 
H  denote  a  separable  Hilbert  space.  Then  a  contraction  A  :  H  — ►  H  is  a 
linear  bounded  operator  such  that  ||T||  <  1.  Now  given  two  contractions  T 
and  A  on  H,  we  set 

Pt(A)  :=  inf { || A/AA/-1 1|  :  M  is  invertible  and  MT  =  TM}. 

The  quantity  pj(A)  is  called  the  T-spectral  radius  of  A. 

Now  given  m  e  H°°  an  all-pass  (inner)  rational  function,  set  H(m)  := 

H2  0  mH2,  and  H  :=  H(m)  0 ...  0  H(m)  (direct  sum  of  p  copies  of  H(rn)). 

For  S  the  unilateral  shift  on  H2 ,  and  for  PH(m)  :  H 2  — ►  H{m)  orthogonal 
projection,  let  S{m)  :=  PH(m^S\H{m)  denote  the  compressed  shift.  We  set 
T  :=  5(m)  ©  ...  0  S(m),  and  U  :=  S  0  . . .  0  5  on  K  :=  H2  0  . . .  ©  H2 
(direct  sums  of  p  copies).  We  can  state  the  following  spectral  version  [19]  of 
the  commutant  lifting  theorem: 

Theorem  8  Taking  T,  U ,  H,  K  as  just  defined,  let  A  :  H  —»  H  be  any 
contraction  such  that  AT  =  T A.  Then  we  have  that 

pr{A)  := 

inf{sr(Z?)  :  B  is  a  commuting  dilation  of  A,  i.e.,  BU  =  UB,  PjjJ9|H  =  APft), 
where  Pn  :  K  — »  H  denotes  orthogonal  projection. 


Remarks.  First  of  all,  note  that  Theorem  7  may  be  derived  as  a  corollary 
of  Theorem  8  by  taking 


m(z) :=  [I 

i=i 


z  —  Zi 

1  -  Z{Z 
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Next  in  [19],  it  is  shown  that  we  can  always  infimize  over  rational  B.  More¬ 
over,  under  the  hypotheses  of  Theorem  8  (in  which  we  are  dealing  with  finite 
block  diagonal  matrices),  one  can  use  gradient  search  procedures  to  compute 
“scalings”  M  arbitrarily  close  to  the  optimal.  Then  if  one  inveitibly  dilates 
such  a  (nearly)  optimal  M  to  an  analytic  matrix  M(z),  one  is  required  to 
check  the  validity  of  matrix  inequality  (16)  for  Mi  :=  M(z{)  (i  =  1  ,...,n), 
i.e.,  via  these  scalings  we  are  reduced  to  a  classical  matricial  Nevanlinna-Pick 
problem. 

Actually  we  have  now  even  found  a  way  to  apply  such  methods  to  the 
structured  singular  value  (in  which  considers  much  more  general  struc¬ 
tured  multiplicative  perturbation  models  as  in  the  work  of  Doyle  and  Sa¬ 
fonov),  which  may  lead  to  an  analytic  procedure  for  performing  /x-synthesis. 
This  will  be  described  in  our  yearly  report  for  AFOSR-90-0024. 

5  Results  in  Interpolation  Theory 

We  have  also  made  some  strides  for  understanding  the  minimal  entropy  so¬ 
lutions  for  a  number  of  interpolation  problems  [15],  [20].  Indeed,  in  studying 
the  spectral  properties  of  the  four  block  operator  and  its  relation  to  classical 
interpolation  theory,  we  have  been  led  to  a  strong  version  of  the  classical 
Parrott’s  theorem.  Parrott’s  theorem  is  one  of  the  key  matrix  extension 
results  and  has  found  numerous  uses  in  control  theory  as  well  as  signal  pro¬ 
cessing.  Our  strengthened  version  in  a  certain  sense  picks  out  an  extension 
which  is  “opposite”  to  the  famous  maximal  entropy  or  central  solution  to 
such  extension  problems.  This  solution  has  a  natural  physical  interpreta¬ 
tion  in  the  waves  through  multi-layered  media  context,  and  we  are  thus 
quite  excited  about  some  of  the  applied  implications  of  our  result.  From  the 
more  theoretical  side,  this  result  can  be  used  in  generalizing  some  beautiful 
results  of  Adamjan-Arov-Krein  on  the  connection  of  the  singular  values  of 
the  Hankel  operator  to  optimal  interpolation  by  functions  with  a  presribed 
number  of  poles  on  the  unit  disc  to  the  four  block  operator  of  H°°  control. 

Further,  based  on  the  strong  Parrott  theorem  we  have  proven  a  strong 
version  of  the  commutant  lifting  theorem  [20]  which  leads  to  an  explicit 
parametrization  of  minimal  entropy  solutions  in  dilation  theory.  This  has 
already  been  applied  to  the  classical  Nevanlinna-Pick  and  Caratheodory 
interpolation  problems. 


6  Other  Results 


During  the  course  of  AF0SR-88-0020,  we  have  have  also  obtained  some  new 
results  on  the  asymptotic  stability  of  time- varying  systems  [16].  Some  work 
has  also  been  done  by  the  Tannenbaum’s  Ph.D.  student,  Hitay  Ozbay,  on 
Ll  optimization  [41].  (Ozbay  showed  that  the  problem  of  minimizing  the 
effect  of  the  worst  disturbance  on  the  output  where  the  signals  have  finite 
one-norm  is  equivalent  to  the  problem  with  bounded  disturbances.) 

Finally,  we  returned  to  the  topic  of  the  studying  the  properties  of  global 
families  of  systems  [23].  This  is  an  area  which  we  plan  to  consider  in  some 
depth  in  AF0SR-90-0024.  Indeed,  previously  our  approach  was  quite  al¬ 
gebraic  and  geometric  (as  exhibited  in  [23]).  We  now  plan  to  use  some  of 
the  analytical  operator  theoretic  techniques  described  above  in  studying  this 
important  problem. 

7  Papers  of  Allen  Tannenbaum  and  Collabora¬ 
tors  under  AFOSR-88-0020 

(All  the  authors  are  listed  in  alphabetical  order  on  the  following  papers.) 

1.  “On  the  Nehari  problem  for  a  certain  class  of  L°°-functions  appearing 
in  control,  II”  (with  C.  Foias),  Journal  of  Functional  Analysis  81  (1988), 
207-218. 

2.  “Weighted  sensitivity  minimization:  general  plants  in  H°°  and  rational 
weights”  (with  T.  Lypchuk  and  M.  Smith),  Linear  Algebra  and  its  Applica¬ 
tions  109  (1988),  71-90. 

3.  “Optimal  sensitivity  theory  for  multivariate  distributed  plants”  (with  C. 
Foias),  International  Journal  of  Control  47  (1988),  985-992. 

4.  “Generalized  interpolation  theory  in  control”  (with  B.  Francis),  Mathe¬ 
matical  Intelligencer  10  (1988),  48-53. 

5.  “On  skew  Toeplitz  operators,  I”  (with  H.  Bercovici  and  C.  Foias),  Oper¬ 
ator  Theory:  Advances  and  Applications  29  (1988),  21-45. 

6.  “A  Poincard-Dulac  approach  to  a  nonlinear  Beurling-Lax-Halmos  theo¬ 
rem”  (with  J.  Ball,  C.  Foias,  and  J.  W.  Helton),  Journal  of  Math.  Analysis 
and  Applications  139  (1989),  496-514. 


25 


n 


7.  “Some  explicit  formulae  for  the  singular  values  of  certain  Hankel  oper¬ 
ators  with  factorizable  symbol”  (with  C.  Foias  and  G.  Zames),  SIAM  J. 
Mathematical  Analysis  19  (1988),  1081-1091. 

8.  “On  the  four  block  problem,  I,”  (with  C.  Foias),  Operator  Theory:  Ad¬ 
vances  and  Applications  32  (1988),  93-112. 

9.  “Iterated  commutant  lifting  for  systems  with  rational  symbol”  (with  C. 
Foias),  Operator  Theory:  Advances  and  Applications  41  (1989),  255-277. 

10.  “On  the  parametrization  of  the  suboptimal  solutions  in  generalized  in¬ 
terpolation”  (with  C.  Foias),  Linear  Algebra  and  its  Applications  124  (1989), 
145-164. 

11.  “On  the  four  block  problem,  II:  the  singular  system”  (with  C.  Foias), 
Integral  Equations  and  Operator  Theory  11  (1988),  726-767. 

12.  “Weighted  sensitivity  optimization  for  nonlinear  systems”  (with  C. 
Foias),  SIAM  Journal  on  Control  and  Optimization  27  (1989),  842-860. 

13.  “Some  remarks  on  optimal  interpolation”  (with  C.  Foias),  Systems  and 
Control  Letters  11  (1988),  259-264. 

14.  “Skew  Toeplitz  approach  to  the  //“''-control  of  multivariate  distributed 
systems”  (with  H.  Ozbay),  to  appear  in  SIAM  J.  Control  and  Optimization. 

15.  “A  strong  Parrott  theorem”  (with  C.  Foias),  Proceedings  of  the  American 
Mathematical  Society  106  (1989),  777-784. 

16.  “Control  of  slowly  time- varying  linear  systems”  (with  E.  Kamen  and  P. 
Khargonekar),  IEEE  Transactions  Automatic  Control  AC-34  (1989),  1283- 
1285. 

17.  “Four  block  problem  :  stable  weights  and  rational  weightings”  (with  P. 
Khargonekar  and  H.  Ozbay),  International  Journal  of  Control  50  (1989), 
1013-1023. 

18.  “On  the  structure  of  suboptimal  H°°  controllers  in  the  sensitivity  mini¬ 
mization  problem  for  distributed  stable  plants”  (with  H.  Ozbay),  to  appear 
in  Automatica. 

19.  “A  spectral  commutant  lifting  theorem”  (with  H.  Bercovici  and  C. 
Foias),  to  appear  in  Transactions  of  the  AMS. 

20.  “On  certain  minimal  entropy  extensions  appearing  in  dilation  theory” 
(with  C.  Foias  and  A.  Frazho),  to  appear  in  Linear  Algebra  and  Its  Appli¬ 
cations. 


26 


21.  “On  spectral  tangential  Nevanlinna-Pick  interpolation”  (with  H.  Bercovici 
and  C.  Foias),  to  appear  in  Journal  of  Math.  Analysis  and  Applications. 

22.  “The  invariant  subspaces  of  a  uniform  Jordan  operator”  (with  H.  Bercovici), 
to  appear  in  Journal  Math.  Analysis  and  Applications. 

23.  “Invariant  theory  and  families  of  systems,”  to  appear  in  a  volume  dedi¬ 
cated  to  Rudolf  Kalman  on  his  60-th  birthday. 

24.  “Operator  theoretic  methods  in  robust  control”  (with  C.  Foias),  to  ap¬ 
pear  in  Signal  Processing  and  Control,  IMA  Series  on  Applied  Mathematics, 
Springer- Verlag,  1990. 

25.  “On  the  singular  values  of  the  four  block  operator  and  certain  generalized 
interpolation  problems”  (with  C.  Foias),  in  Analysis  and  Partial  Differential 
Equations,  edited  by  Cora  Sadosky,  Marcel  Dekker,  New  York  (1989),  483- 
493. 

26.  “On  the  synthesis  of  H°°  optimal  controllers  for  infinite  dimensional 
plants”  (with  H.  Ozbay),  to  appear  in  New  Trends  and  Applications  in 
Distributed  Parameter  Control  Systems,  Marcel  Dekker,  New  York  1990. 

27.  “Spectral  Nevanlinna-Pick  interpolation  theory,”  Proceedings  of  IEEE 
Conference  on  Decision  and  Control,  Los  Angeles,  1987,  1635-1638. 

28.  “On  the  spectra  and  invertibility  of  a  certain  class  of  operators  in  con¬ 
trol”  (with  C.  Foias),  Proceedings  of  IEEE  Conference  on  Decision  and 
Control,  Los  Angeles,  1987,  1338-1342. 

29.  “The  four  block  problem  for  distributed  systems”  (with  C.  Foias),  Pro¬ 
ceedings  of  the  IEEE  Conference  on  Decision  and  Control,  Austin,  Texas, 
1988,  993-998. 

30.  “Spectral  radius  interpolation  and  robust  control”  (with  H.  Bercovici 
and  C.  Foias),  Proceedings  of  the  28th  IEEE  Conference  on  Decision  and 
Control,  1989,  Tampa,  Florida,  916-918. 

31.  “Remarks  on  //“-optimization  of  multivariate  distributed  systems,” 
Proceedings  of  the  IEEE  Conference  on  Decision  and  Control,  Austin,  Texas, 
1988,  985-987. 

32.  “Nonlinear  H°°  theory”  (with  C.  Foias),  to  appear  in  MTNS  Proceed¬ 
ings. 

33.  “Standard  problem  for  distributed  systems”  (with  C.  Foias),  to  appear 
in  MTNS  Proceedings. 


27 


34.  “Robust  control  design  for  a  flexible  beam  using  a  distributed  parameter 
H°°  method”  (with  K.  Lenz,  B.  Morton,  H.  Ozbay,  J.  Turi),  Proceedings  of 
the  IEEE  28-th  Conference  on  Decision  and  Control,  Tampa,  Florida,  1989, 
2673-2678.  1989. 

35.  “On  the  nonlinear  mixed  sensitivity  problem”  (with  D.  Enns,  C.  Foias, 
T.  Georgiou,  M.  Jackson,  B.  Schipper),  Proceedings  of  the  IEEE  28-th  Con¬ 
ference  on  Decision  and  Control,  Tampa,  Florida,  1989,  986-990. 

36.  “Interpolation  theory  in  robust  control”  (with  H.  Bercovici  and  C. 
Foias),  to  appear  in  IFAC  Proceedings,  1990. 

37.  “Spectral  variants  of  the  Nevanlinna-Pick  interpolation  problem”  (with 
H.  Bercovici  and  C.  Foias),  MTNS  Proceedings. 

38.  “A  solution  to  the  standard  H°°  problem  for  multivariable  distributed 
systems”  (with  H.  Ozbay),  Proceedings  of  IEEE  28-th  Conference  on  Deci¬ 
sion  and  Control,  Tampa,  Florida,  1989,  1444-1446. 

39.  “On  approximately  optimal  H°°  controllers  for  distributed  systems” 
(with  H.  Ozbay),  Proc.  of  IEEE  28-th  Conference  on  Decision  and  Control, 
Tampa,  Florida,  1989,  1454-1460. 

40.  “Controller  design  for  unstable  distributed  plants,”  to  appear  in  the 
Proceedings  of  the  ACC,  May  1990. 

41.  (This  paper  supported  under  AFOSR-88-0020  was  done  solely  by  Hitay 
Ozbay.) 

“Z1  optimal  control,”  IEEE  Trans.  Automatic  Control  AC-29,  August 
1989. 


I  * 


Monograph 

J.  Doyle,  B.  Francis,  and  A.  Tannenbaum,  Feedback  Control  Theory,  to  be 
published  by  MacMillan. 


Students  of  A.  Tannenbaum  Supported  by  AFOSR-88-0020 

1.  Hitay  Ozbay  (completed  Ph.D.  in  August  1989) 

2.  Brian  Schipper  (half  time  support  for  masters  degree) 

3.  Evgeny  Berzon  (half  time  support  on  AFOSR-88-0020  and  half  time  on 
AFOSR-90-0024;  he  will  go  to  full  time  support  on  AFOSR-90-0024  in  the 
spring  of  1990  for  his  Ph.D.) 


29 


